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ABSTRACT 

The  Korean  Air  Force  has  determined  that  repairables  management  is  one  of  the 
areas  in  which  attention  could  be  expected  to  lead  to  substantial  improvement  in  the 
efficient  management  of  defense  resources  and  in  maintaining  an  adequate  level  of  force 
effectiveness.  This  thesis  reviews  various  inventory  models  for  the  management  of  re- 
pairable items  and  develops  an  inventory  model  for  the  Korean  Air  Force.  It  discusses 
the  characteristics  of  each  model,  and,  identifies  and  explains  the  differences  in  each 
model  with  respect  to  assumptions,  objectives,  constraints,  and  optimization  methods. 
Also  this  research  evaluates  the  proposed  model  using  sample  data  sets. 
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I.     INTRODUCTION 

A.     BACKGROUND 

The  objective  of  the  Korean  Air  Force  supply  system  is  to  ensure  the  desired  level 
of  operational  availability  for  its  population  of  repairable  equipment  and  maintain  suf- 
ficient stock  levels  for  replacement  components  and  repair  parts  to  support  the  mainte- 
nance of  existing  weapon  systems  and  equipment.  Also,  the  Korean  Air  Force  supply- 
system  is  responsible  for  stocking  a  number  of  items  which  are  failure  prone  and  have 
been  designated  as  repairable.  Because  of  engineering  and  economic  considerations,  at- 
tempts are  made  by  Korean  Air  Force  repair  activities  to  restore  these  items  to  service- 
able condition  whenever  they  fail. 

The  repairable  inventory  system  differs  from  that  of  the  consumable  inventory  sys- 
tem in  two  ways.  First,  the  repairable  system  contains  two  distinct  inventories,  one  of 
which  has  items  that  are  in  a  Ready-For-Issue  (RFI)  state  and  another  which  has  items 
that  are  in  a  Non-Ready-For-Issue  (NRFI)  state.  The  first  inventory  contains  items 
which  are  usable  and  the  second  contains  items  that  must  be  repaired  before  they  can 
be  used.  Second,  the  RFI  inventory  is  made  up  of  a  mixture  of  new  items  and  items  that 
have  been  used,  failed,  repaired,  and  are  ready  to  be  used  again.  This  is  one  reason  why 
repairable  inventory  management  is  a  complicated  process.  [Ref.  1:  p.  2] 

The  Korean  Air  Force  Logistics  Command  (KAFLC)  tries  to  determine  an  appro- 
priate stock  level  for  repairables  based  on  a  primitive  mathematical  model.  This  model 
estimates  stock  levels  using  rules  of  thumb  based  on  expert  opinion.  However,  it  is  very 
difficult  to  determine  an  optimal  stock  level  for  a  repairable  without  using  a  more  so- 
phisticated inventory  management  model.  The  result  of  the  lack  of  such  a  model  is  that 
the  Korean  Air  Force  supply  system  tends  to  maintain  relatively  high  stock  levels  for 
repairable  components.  Although,  repairable  items  account  for  only  a  small  percentage 
of  the  quantity  of  items  stocked  throughout  the  Korean  Air  Force,  they  account  for  a 
large  portion  of  dollars  invested  in  inventory.  The  characteristics  of  high  cost  and  high 
essentiality  that  these  items  typically  possess  make  efficient  inventory  control  difficult 
and  extremely  critical. 

Historically,  most  inventory  models  have  been  developed  for  the  private  sector 
where  the  profit  motive  is  important.  Such  models  consider  the  various  average  annual 
variable  costs  associated  with  managing  inventories  and  strive  to  minimize  the  sum  of 


these  costs.  The  typical  relevant  costs  in  most  inventory  models  are  ordering  cost, 
holding  cost  and  stockout  cost.  Certainly,  these  cost  parameters  are  very  important  in 
the  common  concept  of  inventory  management.  [Ref.  2:  p.  1] 

However,  in  military  supply  systems,  it  is  hard  to  estimate  the  costs  associated  with 
stockouts.  In  addition,  the  military  supply  system  is  not  interested  in  profit 
maximization  or  cost  minimization.  Instead,  such  an  organization  is  usually  interested 
in  maximizing  the  ability  of  its  forces  to  respond  to  any  threat.  An  objective  which 
maximizes  some  measure  of  readiness  with  given  resources  is  therefore  appropriate. 

B.  PURPOSE  OF  RESEARCH 

The  effective  accomplishment  of  each  flying  mission  is  dependent  upon  having  suf- 
ficient numbers  of  aircraft  ready  to  fly  and  to  perform  at  their  fullest  capability.  To 
support  this  goal,  the  Korean  Air  Force  emphasizes  an  extensive  system  of  supply  and 
maintenance  capabilities.  The  Korean  Air  Force  also  tries  to  maximize  the  availability 
of  aircraft  by  quickly  identifying  malfunctioning  parts,  removing  them,  and  rapidly  in- 
stalling replacements  which  have  been  positioned  at  the  base  supply  department. 

The  purpose  of  this  thesis  is  the  development  of  an  inventor}'  model  for  the  man- 
agement of  repairable  items  at  the  wholesale  level  by  the  Korean  Air  Force  Logistics 
Command  (KAFLC).  This  will  be  accomplished  by  the  analysis  of  various  repairable 
inventory  management  models  which  have  been  developed  and  applied  to  the  manage- 
ment of  inventories  of  repairable  items  in  the  United  States  military.  And  also,  a  nu- 
merical analysis  will  be  accomplished  based  on  a  proposed  model.  Finally  the  inventory 
system  parameters,  and  management  concept  will  be  studied. 

C.  PREVIEW 

Chapter  II  of  this  thesis  presents  an  overview  of  the  current  Korean  Air  Force  Lo- 
gistics Command  (KAFLC)  repairable  items  management  process.  This  overview  will 
explain  the  organization  of  the  Korean  Air  Force  Logistics  Command,  describe  the 
present  general  methodology  used  for  repairable  items  management  and  the  current 
overall  approach  to  the  problem. 

Chapter  III  gives  a  detailed  description  of  the  specific  mathematical  inventory- 
models  in  use  today  for  the  management  of  repairables.  Some  explanation  of  their 
mathematical  approach  in  solving  the  inventory  problem  is  given. 

In  Chapter  IV,  the  development  of  a  model  for  the  Korean  Air  Force  repairables 
inventory  management  system  is  described.  Model  formulation  and  solution  techniques 
are  provided. 


Chapter  V  will  show  numerical  results  from  a  computer  program  for  comparing  the 
proposed  model  and  describe  major  components  based  on  the  findings  in  Chapters  III 
and  IV. 

Finally.  Chapter  VI  contains  a  summary  and  recommendations  for  future  study 
based  on  the  proposed  repairables  inventory  management  model. 


II.     REPAIRABLES  AND  THE  KOREAN  AIR  FORCE  LOGISTICS 

SYSTEM 

A.  INTRODUCTION 

The  Korean  Air  Force  didn't  use  any  kind  of  repairables  inventory  management 
model  until  the  late  1960's.  In  the  beginning  of  the  1970's,  the  first  inventory  manage- 
ment concept  was  introduced  by  the  U.S  Air  Force.  The  economic  inventory  manage- 
ment model  was  used  during  this  period. 

The  model  adopted  in  the  1970's  has  been  revised  systematically  several  times  since 
its  adoption  by  the  Korean  Air  Force.  However,  the  usefulness  of  this  inventory  model 
has  decreased  as  the  weapon  systems  used  in  the  Korean  Air  Force  have  become  more 
sophisticated,  and  the  overall  size  of  the  Korean  Air  Force  has  increased.  The  annual 
budget  for  repairables  takes  approximately  60-70  percent  of  the  total  annual  stock  fund 
budget  of  the  Korean  Air  Force,  however  the  essentiality  of  high  quality  repairables  in- 
ventory management  was  not  realized.  [Ref.  1:  p.  1] 

In  this  chapter,  the  Korean  Air  Force  Logistics  system  will  be  described.  The 
chapter  specifically  deals  with  repairables  management,  organization,  system  parame- 
ters, and  mathematical  models  related  to  repairable  item  management.  This  chapter  will 
provide  the  basis  for  the  analysis  and  suggested  improvements  to  the  Korean  Air  Force 
repairables  inventory'  management  system  which  will  be  described  later. 

B.  THE  KOREAN  AIR  FORCE  LOGISTICS  COMMAND  ORGANIZATION 

The  Korean  Air  Force  Logistics  Command  (KAFLC)  is  responsible  for  managing 
the  allocation  of  the  logistics  resources  to  support  maximum  weapon  system  operational 
availability.  It  functions  as  an  intermediate  echelon  command  among  the  tactical  units 
of  the  Korean  Air  Force  under  the  policies  and  planning  of  Air  Force  Headquarters. 
Figure  1  shows  the  organization  of  the  KAFLC.  The  following  paragraphs  will  explain 
the  basic  missions  and  responsibilities  of  the  six  major  divisions  of  the  command. 

The  DM  VI  (Directorate  of  Materials  Management)  is  the  focal  point  of  material 
management  for  the  Korean  Air  Force.  It  procures  all  materials  according  to  its  esti- 
mates of  requirements,  distributes  the  material  to  all  of  the  tactical  units  and  the  sup- 
porting units  under  the  policies  of  the  KAFLC.  Thus,  the  DMM  is  the  equivalent  of  a 
wholesale  level  Inventory  Control  Point  (ICP).  The  DMM  manages  approximately 
210,000  items  which  are  allocated  among  the  hundreds  of  item  managers.    To  manage 


Figure   1.      The  Organization  of  the  Korean  Air  Force  Logistics  Command 


effectively,  DMM  relics  heavily  on  the  Electronic  Data  Processing  Center's  (ED PC) 
computer  services  to  access  historical  transaction  data,  compute  demand  forecasts,  cur- 
rent asset  data.  etc. 

The  DST  (Depot  of  Storage  and  Transportation)  is  the  centralized  warehouse  for 
the  Korean  Air  Force  where  all  procured  and  repaired  materials,  including  consumable1-, 
arc  stored.  They  control  the  physical  inventory  of  the  items  needed  by  the  ultimate  user. 
The  DST  also  deals  with  the  disposal  of  salvage,  scrap,  excess  and  obsolete  material. 
The  other  important  role  of  the  DST  is  transporting  all  materials  to  arrive  at  the  tactical 
unit  when  required. 

The  DME  (Depot  of  Maintenance  and  Equipment)  is  the  in-huuse  depot  level 
maintenance  organization  for  the  Korean  Air  Force.  The  primary  task  of  the  DME  is 
the  overhaul  of  aircraft  on  a  scheduled  basis  (preventive  maintenance).  It  performs  the 
maintenance  of  aviation  end  items  which  arc  beyond  the  capabilities  of  the  the  interme- 
diate and  organizational  maintenance  level  facilities.  These  items  would  be  the  major 
end  items  of  an  aircraft  such  as  the  engine,  fuselage,  gearbox  and  the  direct  ground 
support  equipment  for  the  operation  of  the  aircraft. 


The  D.MEC  (Depot  of  Maintenance,  Electronics  and  Communication)  performs  the 
maintenance  of  radar  equipment  and  its  subassemblies,  ground  communication  equip- 
ment and  weather  forecasting  equipment.  The  DMEC  also  performs  preventive  main- 
tenance. 

The  DMA  (Depot  of  Maintenance  and  Ammunition)  performs  the  maintenance  of 
the  precision  measurement  equipment  (PME),  airborne  equipment  and  armament  of  the 
aircraft,  such  as  laser  guidance  system. 

The  EDPC  (Electronic  Data  Processing  Center)  provides  logistics  software  devel- 
opment and  support,  establishes  job  processsing  standardization,  collection  of  data  re- 
lated to  preventive  and  corrective  maintenance,  provides  analysis  and  guidance  for  each 
of  the  other  major  divisions.  The  major  objective  of  the  EDPC  is  to  assist  each  division 
in  performing  their  specialized  logistics  missions  within  the  Korean  Air  Force  Logistics 
Command. 

C.     REPAIRABLE  ITEMS  MANAGEMENT  SYSTEM 
1.     The  General  System  Overview 

In  the  Korean  Air  Force  an  item  of  supply  is  designated  as  a  repairable  if  it  can 
be  repaired  faster  and  less  expensively  than  it  can  be  procured.  Weapon  systems  in- 
stalled in  aircraft  have  become  increasingly  sophisticated  and  complex.  Hence,  many 
weapon  systems  are  made  up  of  a  number  of  subsystems  which  in  turn  are  comprised 
of  several  repairable  modules.  Often  the  complexity  of  these  individual  modules  are  such 
that  the  personnel  and  equipment  are  not  available  at  the  end  user  level  to  repair  failed 
units.  Consequently,  these  modules  are  designated  as  repairables  and  failed  units  are 
returned  to  designated  repair  activities  for  repair.  Therefore,  repairables  management 
has  become  an  essential  part  of  the  Korean  Air  Force  supply  system  in  terms  of  decision 
making.    Figure  2  shows  the  current  Korean  Air  Force  repairables  management  cycle. 

The  Korean  Air  Force  has  two  kinds  of  maintenance  operations  for  repairables, 
the  base  maintenance  (organizational  and  intermediate  level)  and  the  depot  level  main- 
tenance (DME,  DMEC,  DMA).  It  should  be  noted  that  not  all  depot  level  repairables 
are  repaired  at  the  Korean  Air  Force  facilities.  Because  of  the  complexity,  some  repairs 
are  done  by  the  USAF  under  Foreign  Military  Sales  (FMS)  using  commercial  contrac- 
tors. As  depicted  in  Figure  2,  three  kinds  of  repairables  flows  are  shown.  The  first  is  the 
flow  of  carcasses  to  base  level  maintenance.  The  second  is  the  flow  of  carcasses  to  depot 
level  maintenance.  The  third  level  is  the  flow  of  carcasses  from  the  DST  to  the  USAF 
(the  activities  that  coordinate  with  USAFLC  for  repair  and  procurement  follow  FMS 
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Figure  2.      The  Korean  Air  Force  Repairable  Management  Cycle 


procedures).  Note  that  the  third  flow  implies  the  exchange  of  a  carcass  for  a  new  item 
from  the  USAFLC. 

When  a  failure  of  a  system  or  subsystem  occurs,  the  equipment  is  inspected  to 
determine  the  cause  of  the  failure.  After  isolation  of  the  failure  to  a  component  or  major 
assembly,  further  inspection  is  done  to  determine  if  base  level  repair  is  possible.  In  case 
of  the  failure  of  a  major  end  item,  an  RFI  (Ready-for-Issue)  item  from  the  base  supply 
is  used  to  replace  the  failed  item  and  the  aircraft  returns  to  operational  status.  This  in- 
sures the  operational  availability  of  an  aircraft.  If  possible,  the  failed  item  is  repaired 
by  the  Base  Maintenance  Squadron  (organization  or  intermediate  level  of  repair).  After 
such  a  repair,  the  item  is  released  to  the  Base  Supply  Squadron  in  serviceable  condition. 

For  the  determination  of  the  maintenance  level,  the  Korean  Air  Force  has  a 
special  code  assigned  to  each  replaceable  component  of  a  major  end  item.  Thus,  the 
level  of  maintenance  for  each  component  is  predetermined  by  this  code.  The  code  is  as 
follows.  [Ref.  3:  p.  7] 

First  code  position 

X  expendable  item 

N  repairable  item 

Second  code  position 

D  depot  level  maintenance  item 

F  field  level  maintenance  item  (intermediate) 

B  base  level  maintenance  item 

Third  code  position 

1  depot  level  maintenance 

2  intermediate  level  maintenance 

3  no  maintenance  action 

After  the  Base  Supply  Squadron  has  issued  a  serviceable  item  from  its  stocks 
and  when  base  level  maintenance  has  been  determined  to  be  impossible,  the  Base  Supply- 
Squadron  sends  the  failed  item  to  the  DST.  Of  course,  all  these  actions  are  reported  to 
and  directed  centrally  by  the  DMM.  The  failed  items  from  each  base  are  placed  in  the 
DST  facility  and  await  repair.  These  carcasses  are  turned  over  to  the  depot  maintenance 
facility  (DME,  DMEC,  DMA)  in  batches  under  the  approval  of  the  DMM.    In  other 


words,  all  carcasses  from  each  base  are  stocked  in  the  DST  facility  during  the  quarter. 
At  the  end  of  each  quarter,  all  carcasses  are  sent  to  the  depot  maintenance  facility 
(DME,  DMEC.  DMA).  However,  at  this  time,  item  managers  of  the  DMM  consider 
the  capacity  limit  of  the  depot  maintenance  facility,  l 

2.     Operating  Characteristics  of  Repairables  Management 

As  mentioned  earlier,  there  are  several  depot  maintenance  organizations  in  the 
Korean  Air  Force.  However,  their  concept  of  operation  and  relations  with  the  DMM 
are  identical  for  the  DME.  DMEC  and  DMA.  So,  only  the  operation  of  DME  will  be 
covered  here.  The  maintenance  organization  of  the  Korean  Air  Force  is  restricted  in 
terms  of  its  capacity  to  process  incoming  repairs  and  in  terms  of  the  level  of  technology 
to  deal  with  the  repair  of  complex  systems.  Thus,  the  MRS  (Material  Repair  Schedule) 
and  MRRL  (Material  Repair  Return  List)  are  established  to  manage  these  restrictions. 

a.  Material  Repair  Schedule  (MRS) 

The  MRS  is  the  repair  schedule  for  the  depot  maintenance  facility.  It  de- 
scribes the  type  and  quantity  of  repairable  material  which  can  be  repaired  during  the 
next  fiscal  year  by  repair  depots  operated  by  the  Korean  Air  Force.  In  September  of  the 
fiscal  year,  all  item  managers  estimate  the  demand  for  each  repairable  item  for  the  next 
year,  and  the  total  repair  quantity,  based  on  data  from  the  last  several  years.  The  inte- 
grated results  for  each  item  are  provided  to  the  DME,  DMEC,  and  DMA.  Then,  the 
DME,  DMEC  and  DMA  set  up  their  repair  schedules  after  discussion  and  coordination 
with  DMM.  The  forecasted  maintenance  requirements  which  are  not  covered  by  the 
MRS  are  turned  over  to  the  MRRL.  The  repair  quantity  for  the  MRS  is  set  up  on  a 
quarterly  basis  and  carcasses  are  issued  from  the  DST  based  on  capacity  of  the  DME. 

b.  Material  Repair  Return  List  (MRRL) 

This  is  the  list  of  repairables  for  which  the  Korean  Air  Force  does  not  have 
depot  level  maintenance  capability,  either  due  to  technology  or  capacity  restrictions.  In 
the  case  of  an  MRRL  repair,  the  DST  sends  the  carcasses  to  a  L'SAF  depot  facility. 
Upon  the  receipt  of  the  carcass  from  the  Korean  Air  Force,  the  USAF  returns  a  serv- 
iceable unit  to  the  DST.  Serviceable  units  from  both  the  L'SAF  and  Korean  Air  Force 
depot  facility  are  integrated  at  the  DST  to  make  up  the  serviceable  stocks  which  are 
available  to  support  for  the  tactical  unit. 


1    This  capacity  limit  established  based  on  the  MRS,  see  paragraph  2.a.  above. 


c.     The  Measure  of  Performance 

The  Korean  Air  Force  supply  system  (wholesale  level)  uses  two  measures 
of  performance.  The  first  one  is  fill  rate.  It  is  also  used  by  the  base  level  supply  system. 
The  second  measure  is  supply  response  time.  It  measures  the  length  of  time  elasped  for 
base  backorders  to  be  satisfied  by  the  stock  from  the  DST.  In  other  words,  it  is  the 
length  of  time  from  the  placement  of  an  order  with  the  DMM  to  the  receipt  of  the  or- 
dered item.  Notice  that  the  DMM  uses  this  more  as  a  priority  rule  for  issuing  the  stock 
than  as  a  general  measure  of  performance  for  management  of  the  repairable  items. 
Currently,  the  Korean  Air  Force  assumes  that  any  backorders  for  spares  result  in  an 
aircraft  being  "not  operationally  ready  because  of  supply"  (NORS).  Thus,  the  Korean 
Air  Force  has  developed  codes  which  indicate  the  NORS  condition.  Each  code  provides 
the  maximum  supply  response  time  requirement  to  fill  backorders.  Table  1  summarizes 
the  supply  response  time  requirements  for  each  demand  with  priority  and  Table  2  de- 
scribes the  meaning  of  each  NORS  code.  [Ref.  4:  p.  81] 

Table  1.     SUPPLY  RESPONSE  TIME  REQUIREMENTS 


Priority 

Type  of  order 

Response  time  (day) 

Circumstances 

03 

Express 

3 

G,  K  NORS 

06 

Semi-express 

14 

A,  F  NORS 

13 

Routine 

30 

Routine 

Table  2.     DESCRIPTION  OF  NORS  CODE 


NORS  Code 

Description 

G 

Aircraft  is  in  a  Totally  inoperational  condition 

K 

Radar  malfunctioning 

F 

Operational,  incapable  of  mission 

A 

NORS  is  anticipated 
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D.     MATHEMATICAL  BACKGROUND  OF  REPAIRABLES  MANAGEMENT 

1.     DDR  (Daily  Demand  Rate) 

This  is  the  average  number  of  demands  per  day  and  it  is  computed  following 
formula :- 


365 


365 


where: 

D,  =  demand  on  the  ith  day  of  the  year 

2.  DRP  (Depot  Repair  Percent) 

It  is  computed  in  the  following  manner  for  each  repairable  item  from  the  last 
three  years  of  historical  data: 

DRP   -      RTS    +    XRt's   +    COSD    *  10° 

where: 

RTS  =  repair  this  station 
NRTS  =  non-repair  this  station 
COND  =  condemned 

RTS  refers  to  the  number  of  units  repaired  at  the  Korean  Air  Force  depot. 
Conversely  NRTS  refers  to  the  number  of  units  repaired  by  the  commercial  contractors 
or  by  the  USAF  under  FMS  policies.  Thus,  NDRP  (Non-Depot  Repair  Percentage)  is 
the  complement  of  DRP. 

3.  RCT  (Repair  Cycle  Time) 

The  Korean  Air  Force  has  two  kinds  of  repair  cycle  time.  The  one  is  the  repair 
cycle  time  associated  with  base  maintenance  and  the  other  is  associated  with  the  depot 
level  maintenance.  Usually  RCT  refers  to  the  time  allowed  for  the  depot  level  mainte- 
nance only.  This  time  allowance  is  the  policy  of  the  Korean  Air  Force.  So,  RCT  is 
constrained3  to  be  not  less  than  30  davs  and  not  more  than  120  days. 


2  This  is  computed  by  the  end  of  year.   The  time  period  of  this  computation  starts  from  the 
first  day  of  the  year  to  the  end  of  the  year. 

3  This  is  the  policy  of  the  Korean  Air  Force,  not  a  value  which  is  computed  from  actual 
observations  which  are  recorded  in  the  inventory  management  database. 
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4.  RCQ  (Repair  Cycle  Quantity) 

This  refers  to  the  quantity  of  the  repairable  item  which  must  be  stocked  to  meet 
demands  during  the  repair  cycle  time.  RCQ  is  applied  to  MRS  items  only  and  computed 
by  the  following  formula: 

RCQ    =    DDR    x    DRP    x    RCT 

where: 

DDR  =  daily  demand  rate 
DRP  =  depot  repair  percentage 
RCT  =  repair  cycle  time  (days) 

5.  OLQ  (Operational  Level  Quantity) 

This  is  the  stockage  quantity  which  is  required  to  cover  forecasted  demand 
during  the  time  an  MRRL  item  is  in  the  hands  of  the  USAF  for  repair.  OLQ  is  related 
to  MRRL  as  RCQ  is  to  an  MRS  item.   The  formula  for  OLQ  is  as  follow: 

OLQ   =    DDR    x    60  day 

where: 

DDR  =  daily  demand  rate 

6.  OST  (Order  and  Shipping  Time) 

OST  is  defined  as  the  time  elasped  from  the  initiation  of  a  procurement  or  re- 
pair order,  to  the  receipt  of  the  order.  In  the  case  of  a  procurement,  the  Korean  Air 
Force  constrains  the  OST  to  be  not  less  than  120  days  and  not  greater  than  365  days. 
For  MRRL  repairables,  OST  is  constrained  to  be  not  less  than  220  days  and  not  greater 
than  465  days.  In  both  cases,  the  upper  and  lower  intervals  are  adopted  to  avoid  an 
inventory  stockout.  The  100  days  increment  in  the  OST  is  due  to  the  additional  trans- 
portation time  for  MRRL  items  from  Korea  to  the  continental  U.S.A. 

7.  OSTQ  (Order  and  Shipping  Time  Quantity) 

OSTQ  refers  to  the  quantity  of  repairable  items  needed  to  meet  the  demand 
during  the  order  and  shipping  time.   The  formula  for  OSTQ  is  as  follows: 

for  MRS  items,       OSTQ  =  DDR  x  NDRP  x  OST 
for  MRRL  items,       OSTQ  =  DDR  x  OST 

where: 
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DDR  =  daily  demand  rate 

XDRP  =  non-depot  repair  percentage 

OST  =  order  and  shipping  time 

8.     SLQ  (Safety  Level  Quantity) 

This  refers  to  the  quantity  of  repairables  stocked  to  meet  demand  during  delays 
in  shipping,  delays  in  maintenance,  or  unexpected  increases  in  demand.  The  formula  for 
SLQ  is  as  follows: 


for  MRS  items,        SLQ    =    v3    x    {RCQ    +    OSTQ) 
for  MRRL  items,        SLQ    =    Jl    x    OSTQ 

where: 

RCQ  =  repair  cycle  quantity 

OSTQ  =  order  and  shipping  time  quantity 

9.     RO  (Requisition  Objective) 

Even  with  the  assistance  of  the  existing  computer  system,  it  is  hard  for  the  item 
managers  to  manage  the  more  than  2,000  items  (per  item  manager)  that  are  assigned  to 
them.  Thus,  the  DMM  designates  items  that  have  two  or  more  requisitions  per  year  as 
requisition  objective  items.   The  item  managers  will  monitor  these  items  more  closely. 

Currently,  the  DMM  applies  a  periodic  inventor}'  review  model  which  is  based 
on  a  policy  of  reviewing  and  repairing  at  fixed  regular  intervals  to  bring  inventory  levels 
back  up  to  the  Requisition  Objective  (RO).  Repair  orders  are  placed  at  predetermined 
intervals-*  (the  expected  demand  between  intervals  plus  some  allowance  for  the  variabil- 
ity of  demand  must  be  considered). 

The  requisitioning  objective  consists  of  three  terms.  These  are  the  order  and 
shipping  time  quantity  (OSTQ),  the  repair  cycle  quantity  (RCQ)  and  safety  level  quan- 
tity (SLQ).  The  requisitioning  objective  is  the  sum  of  these  components.  In  Figure  3. 
the  repair  order  quantity  is  simply  the  difference  between  RO  and  the  on  hand  inventory 
at  time  of  review. 

Since  the  system  parameters  such  as  OST,  RCT,  DDR,  DRP  and  NDRP  differ 
at  each  review  period,  the  level  of  the  RO  may  also  be  different  at  the  time  of  each  re- 


•*  This  predetermined  interval  is  the  repair  induction  interval  on  the  MRS.  The  current  length 
of  the  interval  is  2  months. 


13 


RO 
OL 

>» 

f 

/ 

s. 

J 

> 

s                  

s 

) 

N 

S  L 

n  +  2 

n-M 

n  +  G 

v 

Month 

Figure  3.      The  Korean  Air  Force  Rcpairables  Inventory  Periodic  Review 

view.   Inventory  reviews  are  usually  done  every  2  months.  The  formula  for  RO  is  shown 
below: 

for  MRS  items,        RO    =    SLO    +    RCQ    +    OSTQ 


where: 
SLQ    = 
RCQ    = 
OSTQ    =    DDR  x  NDRP  x  OST 


V3    x    (RCQ    +    OSTQ) 
DDR  x  DRP  x  RCT 


for  MRR.L  items,        RO    =    SLO    +    OLO    +    OSTQ 


where: 

SLQ    = 

73    x 

OSTQ 

OLQ    = 

DDR 

x  60  days 

OSTQ    = 

=    DDR 

x  OST 

14 


10.     History  of  RO  Computation 

This  section  shows  the  history  of  the  RO  computation  formula.  The  Korean 
Air  Force  tries  to  determine  optimal  stock  level  for  the  repairables  inventory  manage- 
ment system.  However,  the  results  of  changing  the  RO  formulation  without  adequate 
consideration  of  the  system  parameters  may  result  in  excess  stock  for  some  items  or 
stockouts  for  some  items. 

a.  Before  1973 

RO    =    OLQ    +    SLQ    +    OSTQ 

where: 

OLQ  =  DDR  x  75 
SLQ  =  DDR  x  15 
OSTQ    =    DDR  x  30 

b.  Between  1973  and  1976 

For  MRS.      RO    =    DDR    x    120 

For  MRRL,      RO   =    RCQ   +    SLQ   +    OSTQ 

where: 

RCQ    =    DDR    x    15 


SLQ  =  N  3    x    (RCQ    +    OSTQ) 
OSTQ    =    DDR    x    84 


e.     Between  19~6  and  1979 

RO    =    RCQ    +    SLQ    +    NCQ    +    OSTQ 

where: 

RCQ    =    DDR    x    DRP    x    RCT 
SLQ  =  3(RCQ  +  OSTQ  +  NCQ) 
NCQ    =    DDR    x    NDRP    x    NCT 
OSTQ    =    DDR    x    NDRP    x    OST 
NCT  =  NRTS     Condemn  Time5 


5   This  (NCT)  is  a  time  spent  for  decisioning  a  condemnation  among  NRTS  items. 
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d.  Between  1979  and  1984 

RO    =    RCQ    +    SLQ    +    OSTQ 

where: 

RCQ    =    DDR    x    DRP    x    RCT 

SLQ    =    (DDR    x    NDRP    x    POST)    +    RCQ 

OSTQ    =    DDR    x    NDRP    x    ROST 

e.  After  1984 

For  MRS,      RO   =    RCQ   +    SLQ   +    OSTQ 

where: 

RCQ    =    DDR    x    DRP    x    RCT 


SLQ  =  v'3    x    (RCQ    +    OSTQ) 
OSTQ    =    DDR    x    NDRP    x    OST 


For  MRRL,       RO    =    OLQ    +    SLQ    +    OSTQ 


where: 

OLQ    =    DDR    x    60 


SLQ  =  v'3  OSTQ 

OSTQ    =    DDR    x    OST 

E.     SUMMARY  AND  PROBLEM  STATEMENT 

As  mentioned  earlier,  military  inventor}'  managers  tend  to  want  to  increase  their 
operational  availability  rather  than  minimize  inventory  costs.  However,  in  the  case  of 
the  Korean  Air  Force  repairables  management  process,  the  mathematical  model  is  in- 
sufficient to  adequately  model  the  stochastic  characteristics  of  the  inventory  manage- 
ment problem.  In  the  RO  (Requisition  Objective)  computation  it  does  not  specify  the 
wear-out  rate  or  the  regeneration  rate,  which  are  important  parameters  in  dealing  with 
repairables  inventory  management. 

The  DDR  (Daily  Demand  Rate)  computation  is  the  mean  of  a  stochastic  process. 
However,  the  current  inventory  model  treats  demand  as  if  it  were  deterministic.  In  the 
formula  for  SLQ  (Safety  Stock  Quantity),  they  considered  the  stochastic  nature  of  de- 
mand but  the  actual  probability  distribution  is  not  considered  at  all.  As  evidence  of  the 
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model's  shortcomings,  the  Korean  Air  Force  repairables  management  tends  to  stock 
large  quantities  of  spare  parts  and  components  in  order  to  avoid  stockout  situations. 


17 


III.     OVERVIEW  OF  REPAIRABLE  MANAGEMENT  MODEL 

A.  INTRODUCTION 

Inventory  systems  such  as  the  one  described  in  the  previous  chapter  can  be  classified 
rather  broadly  as  a  multiechelon  system  with  repair.  Analytical  studies  of  multiechelon 
systems  have  shown  the  computation  in  such  models  to  be  hard,  so  either  simplifying 
assumptions  or  approximations  are  necessary.  The  introduction  of  the  repair  aspect  into 
the  system  certainly  complicates  the  structure  of  a  model  even  at  single  echelon  levels. 

Suppose  a  repairable  system,  consisting  of  one  Inventory  Control  Point  (ICP),  one 
stock  point  and  one  overhaul  and  repair  activity  controls  the  inventory  of  a  single  item. 
Demands  from  various  tactical  units  are  placed  only  at  the  stock  point.  The  system  has 
continuous  updating  of  records,  i.e.,  transaction  reporting.  When  items  wear  out  or  fail, 
the  tactical  unit  can  either  scrap  the  item  or  return  it  to  the  depot  repair  facility  (DME, 
DMA,  DMEC).  After  inspection,  the  depot  can  either  scrap  the  item  or  repair  and  re- 
turn it  to  the  RFI  inventory. 

In  this  chapter,  we  will  discuss  various  repairable  inventory  models  based  on  two 
distinguishing  characteristics.  One  is  a  deterministic  model,  the  other  is  a  stochastic 
model.  In  case  of  deterministic  model,  we  will  describe  two  models,  "continuous  sup- 
plement" model  and  "substitution"  model.  For  the  stochastic  case,  we  will  describe  five 
models  based  on  current  usage  in  U.S  military.  Each  description  will  include  model 
formulation  and  solution  techniques. 

B.  DETERMINISTIC  INVENTORY  MODEL 

These  deterministic  inventory  models  were  developed  by  David  A.  Schrady  and 
others  between  1966  and  1967.  These  deterministic  inventory  models  are  described  in  a 
research  report  for  the  U.S  Navy  repairable  inventory  management  system  [Ref.  5,  6, 
7].  The  first  model  calls  for  the  repair  activity  to  induct  a  batch  of  carcasses  when  the 
XRFI  inventory  level  reaches  a  certain  point.  With  a  deterministic  system,  this  rule  in- 
sures regularly  spaced  repair  inductions  of  fixed  batch  size,  which  is  a  simple  process  for 
the  repair  activity  to  manage.  Repair  batch  sizes  can  be  maintained  if  there  is  contin- 
uous supplementing  of  repaired  RFI  items  with  new  procurement.  The  procurement 
order  will  be  issued  with  enough  lead  time  so  that  when  the  on-hand  RFI  stock  drops 
to  zero,  a  procurement  quantity  will  arrive.  The  two  inventories,  RFI  and  NRFI,  will 
have  time  histories  as  shown  in  Figure  4.  Note  that  in  this,  the  "continuous  supplement" 
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model,  the  procurement  trigger  is  in  the  RFI  inventory  and  the  repair  trigger  is  in  the 
NRFI  inventory.  [Ref.  6,  7:  pp.  9-10,  pp.  392-393J 

The  second  model  is  suggested  by  noticing  that  there  is  a  trade-ofT  between  stock 
held  in  RFI  condition  and  stock  held  in  NRFI  condition.  NRFI  carcasses  are  a  legiti- 
mate resource  because  they  require  only  repair  to  be  restored  to  full  usefulness.  But  the 
value  of  this  resource  is  less  than  the  value  of  the  RFI  resource  by  at  least  the  cost  of 
labor  and  replacement  parts.  Thus,  if  inventory  is  to  be  held  in  the  system,  it  would  be 
slightly  cheaper  to  hold  it  in  NRFI  condition  than  in  RFI  condition.  [Ref.  6,  7:  p.  10, 
p.  393] 

Notice  that  the  "substitution"  policy  minimizes  the  RFI  inventory  while  the  "con- 
tinuous supplement"  policy  minimizes  the  NRFI  inventory  [Ref.  6:  p.  13].  The  next 
section  describes  the  mathematics  of  these  two  models. 

1.     Analysis  of  The  Continuous  Supplement  Policy  Model 
a.     dotation 

Under  the  assumption  of  deterministic  demand  and  lead  times,  it  is  not 
necessary  to  maintain  a  safety  stock.  Whenever  the  RFI  inventory  reaches  zero,  a  pro- 
curement quantity  will  arrive.  Thus,  the  system  is  never  out  of  stock.  Between  pro- 
curement arrivals,  depleting  RFI  stock  is  replenished  with  repaired  items.  Also,  a 
repetitive  system  will  be  established  regardless  of  the  initial  provisioning  policy,  so  that 
it  is  sufficient  to  analyze  only  one  cycle  to  determine  system  characteristics.  Figure  4, 
shows  basic  characteristics  of  the  "continuous  supplement"  policy  model.  The  basic 
notation  as  follows  [Ref.  5:  pp.  12-22]: 

D  =  annual  demand  rate 

Tx  =  procurement  cycle 

7~:  =  repair  cycle 

r,  =  procurement  lead  time 

fj  =  repair  lead  time 

r0  =  field  recovery  rate 

r2  =  overhaul  and  repair  recover}'  rate 

Q.  =  procured  lot  size 

Q2  =  repaired  lot  size 

hx  =  holding  cost  for  RFI 

h2  =  holding  cost  for  NRFI 

X-  =  reorder  point  for  procurement 
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X2  =  reorder  point  for  repair 

n  =  number  of  repair  cycles  per  procurement  cycle,      n  =  7,  /  T2 

Ax  =  fixed  procurement  cost  per  one  cycle 

A2  =  fixed  repair  cost  per  induction 

b.     Model  Formulation 

The  objective  of  the  model  is  to  determine  an  optimal  procurement  quan- 
tity Q\  and  an  optimal  repair  quantity  Q\  that  minimize  the  average  annual  cost.  These 
optimal  values,  coupled  with  respective  reorder  points,  X\  and  X2,  constitute  an  optimal 
operating  doctrine.  To  determine  the  average  annual  variable  cost,  the  costs  per  repair 
reorder  cycle  must  be  investigated.  The  product  of  per  cycle  costs  and  the  number  of 
cycles  per  year  will  then  yield  the  average  annual  cost.  The  total  cost  for  any  given  cycle 
r,  is  the  sum  of  procurement,  repair  and  holding  costs. 

Since  there  is  only  one  procurement  per  reorder  cycle,  the  variable  order 
cost  is  the  actual  cost  of  the  items  ordered  and  can  be  expressed  as  QXCV  The  fixed 
procurement  cost  for  one  cycle  is  Ax  as  mentioned.  To  determine  the  repair  costs  per 
cycle  it  is  necessary  to  compute  the  number  of  repair  cycles  per  procurement  cycle  ( 
n  =   7",  /  T2,  n  is  an  integer).    In  Figure  4: 

h  =  'i  +  h 

where: 

h  =  KT2  ,  0<K<\ 

t2  =  {\    -    K)T2 

The  cost  of  items  repaired  per  procurement  cycle  is  C2Q2  and  the  fixed  repair  cost  is 
A2n,  where  A2  is  defined  to  be  the  fixed  repair  cost  per  induction.  The  holding  costs  per 
procurement  cycle  will  be  h{AT  +  h2AT2,  where  AT  is  the  area  under  RFI  curve  and 
AT2  is  the  area  under  the  NRFI  curve.  To  compute  the  area  under  RFI  curve,  consider 
Figure  5  showing  RFI  inventory  for  one  procurement  cycle.  Since  D  is  constant  and 
known,  the  area  of  U,  AUt  is  given  by: 

Av=   y  /,  (2,  +  a) 
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Figure  4.      Repair  Cycle  and  System  for  "Continuous  Supplement"  Model 
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where: 

h    =    KT2 

a   -   Qx   -   Dtx 

The  formula  will  reduce  to: 


Av=  (/,e,  -  Dt])  /  2 


Since  n  is  an  integer,  the  area  of  V  is  the  sum  of  n-1  trapezoids  each  having 
one  side  of  length  T2.  To  determine  the  area  of  each  trapezoid,  a  recursive  relation  is 
developed.    Let: 

Av.  =  area  of  the  ith  trapezoid  in  V  (see  Figure  4) 

=  Av>  +  Av>  +  +^v. 

where: 

AVl  =y  U2Q2  +  2 (ft  -  hD)  -  DT2-\ 

Ay2  =   y  U402  +  2(0,  -  t,D)  -  3Z)72] 
therefore,  the  general  formula  for  the  area  of  V: 

Ay,  =   y  T2\2iQ2  +  2(<2i  -  f,  D)  -  (2/-  1)Z>7"2]       /=  1,2,....,*  -  1. 

n-l 

So  the  area  of  V  is.  Av—  Y.AV\ 

1=1 

^k  =   4-  T2[2Q2^i   +  2(«-l)(&  -  f,D)    -    072(21-  1)1 

=   y  r'£&  » (»  -  1)  +  2(«  -  1)&  -  2(«  -  l)r,Z>  -  w  -  1)  +  Z>r2(/i  -  1)] 

=    jr2[2(«-l)(2,  -  f,D)  +  («2  -  n)Q2  -  (n-lYDT2-] 
Finally,  area  of  W  in  Figure  4,  Aw  is  easily  calculated  as  follow: 


A  iv   —      ->  h 
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where: 

t2  =  (1  -  k)T2 

b  =  n02  +  &  -  t,D)  -  («  -   \)DT2 

Thus,  the  area  of  W: 


*w   =    jT2[nQ2    +    Qx    -    ^Z>    -    (if  -  \)DT{] 


The  total  area  under  RFI  curve  is: 


AT  =  Av  +  Av  +  A 


w 


=  1&  -    -jDt;  +  y  r2{2(/f-l)(&)  -  ZH  +  («2  -  »)&  -  («2  -  n)DT2) 

Since  n  is  an  integer  and  the  build  up  rate  of  XRFI  items,  r0D,  is  constant, 
the  area  under  the  XRFI  curve  is  simply  n  times  the  area  under  the  repair  cycle  curve 
is  T202l2rv  So  the  XRFI  holding  cost  per  procurement  cycle  is  nT2Q2h2\2r2.  The  total 
cost  per  procurement  cycle  becomes: 

A2TX           C20,Tx  h2Q2Tx 

KTl  =  2i  Q  +  A\  +   -J—   +   —j +   — j — Yl  +  hlA^ 

The  total  average  annual  cost  is  then: 

K        ^_  Q^  A±_  A2       C2Q2  h2Q2  h}ATi 

T,  Tx       +     Tx     +     T2+     T2  2/-2      +       r, 

Simplifying  the  total  average  annual  cost: 

A.RD          rDn                             h2Q2           hxQx  hxQ2 

A  =  +  C,^/)  +  ^2A-2r0Z)  +    -y—   +    —2 hxkQ2  +    — y- 

where: 

T2  =  ft  /  r0r2D 

TJ  T2  =  n  ,      n  =   1,2,3,.... 

f,  =  A72  =  *£  /  r0r2Z)     for,    0  <  k  <  1 

f,  +  fa  =  T2 

r,  -  fi,  / 1  -  ^£>  =  QJRD 
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R  =  1  -  r0r2 

Note  that  the  terms  CXRD  and  C2r0r2D  are  independent  of  <2i  and  Q2,  and  hence  do  not 
affect  the  operating  doctrine.  Therefore  it  is  appropriate  to  redefine  the  average  annual 
cost  of  ordering,  repairing  and  holding  as  follow: 

A,RD  A2r0r2D  h202  A,e,  .  ..  h,Q2 

k  =  ~qT  +  ~~QT  +  ~  +T  ~  h'kQl  +  ~ 

From  differential  calculus  the  optimal  Qx  and  Q2  will  be  those  that  satisfy  the  equations, 

SK    _  0        J£.  .  0 


3(2, 

Taking  the  partial  derivatives  results  in 


dK               A,RD           h, 

eg,           Q\         2 

dK 
cQ2 

A2r0r2D             h2                         /t, 

Ql            2>'i                  2 

Solving  the  equation  to  get  the  optimal  procurement  and  repair  quantities, 


.  =         2.4  .RD  .  =      J         2A2r0r20D 


hx        '  ^l        V   h2  +  r2h,{\  -2k) 

2.     The  Substitution  Policy  Model 
a.     Notation 

Notice  that  this  model  minimizes  the  RFI  inventory.  Figure  5  shows  basic 
characteristics  of  the  "substitution  "  policy  model.  The  basic  notation  as  follows  [Ref. 
6,  7:   pp.  13-17,  pp.  393-396]: 

QP  =  procurement  quantity 

QR  =  repair  batch  size 

d  =  demand  rate 

(1  —  r)  =  scrap  rate,  r  is  the  recovery  rate  measured 

tP,  tR  =  procurement  and  repair  lead  times 

AP  =  fixed  procurement  cost  per  order 
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An  =  fixed  repair  batch  induction  cost  per  batch 
ht  =  RTI  holding  cost 
h2  =  NRFI  holding  cost 
T  =  system  cycle  time 

b.     Model  Formulation 

Referring  to  Figure  5,  define  Tc  as  the  time  period  during  which  inductions 
are  suspended  and  the  repair  is  simply  accumulating  NRFI  carcasses.  It  may  be  deter- 
mined that: 

Ta  =  (Qp+  QR)ld 

Next,  let  n  be  the  number  of  inductions  per  cycle.  In  general,  it  will  not  be  possible  to 
insure  that  the  last  induction  before  T0  begins  will  reduce  the  NRFI  stock  to  zero.  The 
residual  NRFI  stock  when  T„  begins  will  be  some  fraction  of  the  net  loss  in  NRFI  items 
per  repair  cycle,  where  the  repair  cycle  is  the  time  between  regular  consecutive  in- 
ductions. The  net  loss  o[  NRFI  inventory  over  the  period  between  successive  in- 
ductions: 

rd{QRjd)  -  QR  =    -QR{\-r) 

Thus,  the  residual  when  Ta  begins  will  be  some  fraction  of  QR{\  —  r),  call  it  fiQR{\  —  r) 
where  0  <  /?  <  1.  While  the  /?  factor  is  unavoidable  due  to  the  requirement  that  there  be 
an  integral  number  of  repair  batches  in  each  cycle,  we  shall  set  /?  =  0  in  subsequent  de- 
velopments. 

Now  the  first  batch  after  inductions  are  resumed  takes  the  amount  QR  from 
the  NRFI  inventory.  Subsequent  inductions  cause  a  net  reduction  of  only  QR{\  —  r) 
items.  Thus,  the  amount  of  NRFI  items  available  for  the  (n-1)  inductions  are 
rdT0  —  QR  or  QR(r  —  1)  +  rOP.    Dividing  the  last  expression  by  QR(\  —  r)  yield: 


QR(r-\)  +   rOP  _^ 

QR(\-r)  +ll-r 


(«-!)=  „  77    ...  =   -  1  +  (-T-t)QpQr 


or 


r         Qp 


Scrop 
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Figure  5.      Repair  Cycles  and  System  for  "Substitution"  Model 
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The  system  cycle  time  is  given  by: 

"Qr  +  Qp  Qp 


T  = 


(1  -  r)d 


where: 

From  the  Figure  7,  the  area  under  the  RFI  curve,  over  one  complete  cycle  Au  is  then: 

A]  =  ~2d{T^7)lQRQp  +  (17J1)^] 

Referring  again  to  Figure  7  the  area  the  XRFI  curve: 

^2  =  y{1^7)LQrQp+  Q2p] 

Now  we  get  total  cost  per  unit  times  expression: 

rr  d(\  —  r)  ARrd  h->r  i  _  ,. 

"T"  =  Ap-h; —  +  -75—  +  -4-C&  +  {±tl')Qp'\  +  hi'LQp  +  Qr\ 

1  SdP  V-R  l 

where: 

Qp  +  "Qr  Qp 


T  = 


d{\-r) 
The  total  cost  under  substitution  model: 

TC  =  AP  +  nAR  +    ^Y^7IQpQr  +    1TJLQp']  +    ^2d\~7lQpQR  +  Q& 

The  optimal  order  quantities  are  obtained  by  setting  the  partial  deriatives, 
with  respect  to  QP  and  QR: 

cTC    =  Q         17_C    =  Q 


6QP  6QR 

Next  we  get  new  formula  for  optimal  order  quantities: 
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cTC 
8QP 

Apd{\-r)           h,{\-r) 

Q\                 2 

h2r 

2 

cTC 
SQr 

ARrd           h\Y            h2r 

Ql      '      2      '      2 

The  optimal  order  quantities  for  repair  and  procurement: 


/     2APd(\-r)  ,  /     2A2d 

Up  ~  a  /  u  i\  _  .a  _l   u  ~  '  ^«  ~ 


Ai(i  -0  +  V  '  V  ^  +  ^2 

3.     Summary  and  Evaluation 

An  important  result  of  deterministic  inventory  model  analysis  is  the  realization 
that  economic  order  quantities  do  not  vary  directly  with  demand  but  as  the  square  root 
of  demand.  This  nonlinear  relationship  could  explain  why  many  organizations  experi- 
ence inventory  problems.  Intuitive  inventory  policies  are  undermined  by  difficult-to- 
grasp  nonlinearities.  [Ref.  8:  p.  141] 

Although  the  models  described  in  this  section  are  evidently  applicable  to  re- 
pairables  management,  they  are  fraught  with  many  limitations.  The  reasonableness  of 
model  assumptions  and  their  sensitivity  to  actual  conditions  determines  the  utility  of  any 
particular  model.  For  example,  the  two  models  assume  that  the  demand  for  a  repairable 
item  is  known  with  certainty.  Also  stockouts  either  are  not  permitted  to  exist  or  are 
backordered  and  satisfied  when  replenishments  are  received. 

In  terms  of  the  Korean  Air  Force  repairable  management  system,  the  models 
described  here  are  difficult  to  adapt  to  the  Korean  Air  Force  system.  This  is  because 
of  demand  variability,  relatively  long  procurement  lead  time,  repair  capacity  and  capa- 
bility, a  number  of  backorder  situations,  and  current  system  management  concept.  The 
other  important  point  is  that  the  Korean  Air  Force  often  encounters  NORS  conditions 
which  results  in  decreased  operational  availability.  However,  these  two  deterministic 
models  do  not  deal  with  stockouts.  These  reasons  didn't  allow  to  use  for  the  Korean 
Air  Force  repairable  management  system. 

C.     STOCHASTIC  INVENTORY  MODELS 

As  the  weapon  systems  installed  in  modern  aircraft  become  increasingly  sophisti- 
cated and  complex,  repairable  items  represents  an  important  subset  of  the  total  inven- 
tory of  items  which  are  managed  by  the  military  supply  system. 
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In  general,  repairable  items  are  supported  by  a  two-echelon  inventor}'  and  repair 
system.  When  a  repairable  item  fails  at  the  base  level,  it  is  returned  to  base  supply  and 
a  new  serviceable  unit  is  issued  from  extra  RFI  stock  of  the  base.  If  possible,  the  failed 
item  is  then  repaired  by  the  base  maintenance  organization  and  then  stored  in  RFI 
conditions  at  base  supply.  Sometimes,  however,  the  failed  item  must  be  returned  to  the 
repair  depot  where  more  sophisticated  equipment  and  specialized  skills  are  available  to 
repair  it. 

When  we  consider  the  condemnation  of  repairables,  there  should  be  an  inflow  of 
newly  manufactured  items  to  depot  level  supply  should  look  at  the  material  flow  of  all 
echelons,  and  then,  place  a  replenishment  order  to  resupply  the  condemned  repairables. 
At  times,  forward  base  locations  will  be  supplied  from  another  closely  located  base.  This 
is  called  lateral  resupply.  For  other  items,  a  manufacturer  may  provide  both  the  source 
of  procurement  for  new  assets  and  the  source  of  repair  for  failed  items. 

The  decisions  concerning  which  maintenance  levels  will  repair  the  failed  items  sub- 
sequently affects  the  supply  support  provided  at  the  organizational  level.  If  an  item  is 
repaired  at  the  organizational  level,  repair  equipment  and  maintenance  personnel  must 
be  made  available  at  that  level.  However,  if  the  item  is  not  repairable  at  that  level,  then 
the  question  is  whether  the  item  can  be  replaced  at  the  organizational  level.  Organiza- 
tional level  repairables  will  normally  be  transferred  to  the  next  higher  echelon  of  repair 
if  the  repairs  cannot  be  accomplished  at  the  organizational  level.  Similarly,  repairables 
which  cannot  be  repaired  by  the  intermediate  level  are  usually  sent  to  the  depot  level  for 
repair. 

The  difference  between  stochastic  models  and  deterministic  models  is  that  if  de- 
mands or  future  requirements  are  uncertain,  then  regardless  of  the  stockage  policy 
adopted  there  is  generally  some  probability  that  available  stock  levels  will  be  insufficient 
to  meet  demand  [Ref.  9:  pp.  187-188].  In  this  section,  we  will  describe  five  stochastic 
models  based  on  current  usage. 
1.     Model  Analysis 

a.     Analysis  of  ME  TRIC  Model 

(1)  Background.  The  METRIC  model  was  developed  over  a  period  of 
years  by  a  research  group  at  the  RAND  corporation  and  presented  in  the  literature  by 
Sherbrooke  (1968)  and  extened  by  Muckstadt  (1973).  METRIC  was  developed  with  the 
ultimate  goal  of  implementation  and  a  slightly  modified  version  of  METRIC  was  actu- 
ally implemented  by  the  USAF.  [Ref.  10.  11,  12] 
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METRIC  is  a  mathematical  model  of  a  base-depot  supply  system  in 
which  item  demand  is  assumed  to  be  compound  Poisson  with  a  mean  value  estimated 
by  a  Bayesian  procedure.  The  model  is  also  capable  of  determining  base  and  depot  stock 
levels  for  stock  levels  of  a  group  of  recoverable  items.  Its  governing  purpose  is  to  opti- 
mize system  performance  for  specified  levels  of  system  investment.  It  is  designed  for 
application  at  the  weapon-system  level,  where  a  particular  line  item  may  be  demanded 
at  several  bases  and  the  bases  are  supported  by  one  central  depot.  The  support  depot 
may  vary  by  item  as  in  the  item-manager  system  or  it  may  be  fixed  as  in  the  weapon- 
system  storage  site  concept.  The  general  purposes  of  this  model  are  optimization,  re- 
distribution and  evaluation.  [Ref.  13:  p.  2] 

The  basic  METRIC  model  considers  a  two-echelon  system  in  which 
independent  bases  (lower  echelon)  are  supported  by  a  repair  depot  (upper  echelon).  At 
the  occurrence  of  a  failure  (it  could  be  more  than  one),  the  failed  item  is  either  replaced 
by  available  base  stock  or  back  ordered  if  the  base  stock  is  not  available.  The  item  is 
inspected  to  determine  the  extent  of  repair  required.  If  the  repair  can  be  made  at  the 
base,  the  unrepaired  item  enters  base  repair.  If  the  item  cannot  be  repaired  at  the  base 
level,  it  is  shipped  to  the  depot.  Upon  shipping  the  item  to  the  depot,  the  base  places 
an  order  with  the  depot  for  a  replacement,  so  that  the  inventory  position  for  item  i  at 
base  j  can  be  maintained. 

(2)     Mathematical  Assumptions 

1.  System  Objective  of  Minimizing  the  Expected  Number  of  Backorders:  The  objective 
will  be  to  minimize  the  sum  of  expected  backorders  on  all  recoverable  items  at  all 
bases  pertinent  to  a  specific  weapon  system.  Thus,  unless  all  bases  are  identical, 
the  expected  number  of  backorders  will  vary  by  base.  Expected  backorders  take  a 
fixed  period  of  time  and  add  together  the  number  of  days  on  which  any  unit  of  any 
item  at  any  base  is  backordered.  In  order  to  minimize  the  expected  number  of 
backorders  one  divides  this  number  by  the  length  of  the  period  and  taking  the  ex- 
pected value.  This  yields  a  number  that  is  independent  of  period  length.  This  is 
the  value  we  seek  to  minimize.  [Ref.  14:  pp.  126-131] 

2.  Compound  Poisson  Demand:  We  assume  that  demand  for  each  item  is  described  by 
a  logarithmic  Poisson  process.  Compound  Poisson  processes  are  generalizations 
of  Poisson  processes;  the  compound  processes  allow  the  flexibility  of  incorporating 
more  parameters,  yet  retain  the  simple  analytical  properties  of  the  Poisson.  The 
logarithmtic  Poisson  is  obtained  by  considering  batches  of  demands  where  the 
number  of  batches  follows  a  Poisson  process  and  the  number  of  demands  per  batch 
has  a  logarithmic  distribution. 

3.  Demand  is  Stationary  over  the  Prediction  Period:  It  is  assumed  that  the  distribution 
of  demand  over  some  future  period  of  interest  is  stationary. 
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4.  Repair  Decision  Deponds  on  the  Complexity:  The  assumption  is  that  the  decision  to 
repair  a  unit  at  base  level  or  send  it  to  the  depot  is  a  function  only  of  the  type  of 
malfunction  and  the  base  maintenance  capability. 

5.  Lateral  Re  supply  is  Ignored:  When  a  unit  is  shipped  from  base  to  depot  for  repair, 
a  serviceable  replacement  will  be  resupplied  from  the  depot  if  possible.  If  the  depot 
has  no  unit  on  the  shelf,  the  base  must  wait  until  a  unit  emerges  from  depot  repair. 

6.  System  is  Conservative:  Consider  a  particular  stock  item  demanded  at  base  j.  We 
assume  that  a  unit  of  stock  has  a  probability  r,  of  being  repairable  at  base  j,  with 
a  repair  time  drawn  at  random  from  a  base  repair  distribution  with  mean  A/,  a 
probability  (1  —  r;)  of  being  depot  repairable,  with  an  order  and  shipping-time 
distribution  having  mean  0,  and  a  depot  repair  distribution  having  mean  D.  This 
implies  that  there  are  no  condemnations  or  that  the  system  is  conservative,  as  the 
name  "recoverable  item"  suggests.  A  higher  condemnation  rate  usually  indicates 
that  the  item  should  be  redesigned.  The  condemnation  rate  must  be  considered  for 
procurement  purposes,  but  the  procurement  process  is  not  considered  in  the 
METRIC  optimization. 

7.  The  Depot  Does  Sot  Batch  Units  of  Recoverable  Item  for  Repair:  The  model  as- 
sumes that  depot  repair  begins  when  the  repairable  base  turn-in  arrives  at  the  de- 
pot. This  appears  to  be  a  reasonable  approximation  to  current  depot  scheduling 
practice.  Since  METRIC  economizes  by  buying  fewer  high-cost  items,  these  are 
the  items  which  are  most  likely  to  be  in  short  supply. 

S.  Demand  Data  from  Different  Bases  Can  Be  Pooled:  We  assume  that  demand  from 
the  several  bases  can  be  pooled  in  some  manner  so  that  a  composite  initial  estimate 
of  demand  per  flying  hour  can  be  obtained.  The  pooled  estimate  can  be  obtained 
by  a  simple  averaging  technique  or  a  more  sophisticated  procedure  such  as  expo- 
nential smoothing.  METRIC  multiplies  this  number  by  the  flying  hours  per  month 
for  each  item. 

r3j  Model  Formulation.  By  mathematical  assumptions,  demand  in  this 
system  follows  a  compound  Poisson  process.  A  compound  Poisson  process  may  be 
thought  of  as  a  series  of  customers  who  arrive  following  a  Poisson  process,  each  of 
whom  can  demand  an  amount  that  is  independently  distributed  according  to  a  com- 
pounding distribution.  Assume  that  item  i  is  stocked  at  each  of  base  j,  and  the  cus- 
tomers who  place  demand  for  the  item  at  each  base  have  a  known  mean  arrival  rate  of 

/,.  j  =  1.2.3 ,J.    When  a  customer  arrives  at  base  to  place  one  or  several  demands,  he 

turns  in  an  equal  number  of  carcasses.  These  carcasses  can  be  repaired  at  base  level  with 
probability  of/-,.,  while  (1  —  rv)  is  the  probability  that  they  must  be  repaired  at  the  de- 
pot. The  arrival  of  carcasses  from  base  j  at  the  depot  is  described  by  a  Poisson  process 
whose  mean  is  (1  —  r/y)  times  the  mean  of  the  Poisson  customer  arrival  process  at  the 
base  j.  Therefore,  the  total  demand  at  the  depot  for  item  i  is  compound  Poisson,  with 
mean  customer  arrival  rate: 
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7=1 

where: 

/.,j  =  mean  arrival  rate  fot  item  i  at  base  j 

Let/,  be  the  mean  demand  per  customer  at  base  j.  Then,  the  mean 
depot  demand  rate  for  unit  i  is: 

j  J 

e  =  EVtiU  -  nj)  =  2>//i  -  r„) 

;=i  ;=1 

where: 

9,,  =  the  mean  rate  for  item  i  at  base  j 

In  the  special  case  of  the  logarithmic  Poisson  process,  the  probability 
that  x  customer  demands  are  in  the  repair,' supply  process  is  negative-binomial  with  pa- 
rameters q  and  K  (Note:  K  =  AT  I  \nq  where  ).  is  mean  customer  arrival  rate  and  T  is 
average  resupply  time). 

P(x  I  XfTg)    =    (K+x-n        iq  -  "  —K 

(K  —   ly.xlq 

x  =  0,1,2, ,        q>  1,       A">0 

where: 

q  =  the  variance  to  mean  ratio  of  P(x  |  XVT^ 

K=  Wvl(q  -   1) 

T{j  =  average  resupply  time  for  a  demand  for  item  i  at  base  j. 

A,j  =  base  repair  cycle  time 

D  =  depot  repair  cycle  time 

Oy  =  the  order  and  shipping  time 

S(S0)  =  average  depot  delay- 
Then: 

Ty  =  ryAy  +  (1   -  ry){Oy  +  S(S0)} 
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Since  it  takes  an  average  of  D  time  units  for  an  arrival  to  complete 
the  repair  process,  the  probability  distribution  of  the  number  of  units  in  the  depot  repair 
cycle  is  compound  Poisson  with  mean  of  AD.  Hence,  the  expected  number  of  units  back 
ordered  at  the  depot  is: 

BO(S0\AD)  =    Yj(x  ~  S0)P(x\AD) 

x>S0 

As  mentioned  in  the  beginning,  the  objective  of  METRIC  is  to  min- 
imize the  sum  of  backorders  for  all  item  i  and  for  all  bases  j  within  a  budget  constraint. 
Thus,  the  METRIC  problem  will  be  represented  as  follows: 

/     j 


minimize  ^  ^BO^S^,  Sy) 


/    j 
subject  to  2,  /    CjSjj  <  C 

StJ  >  0,     1  <  i  <  I,     0<j<J 

where: 

S,;  -  the  decision  variables,  stock  for  item  i  at  base  j 

C  =  the  total  amount  of  budget  available 

C,  =  the  cost  of  item  i 

S,0  =  the  depot  stock  for  item  i 

(4j  Solution  Technique.  The  METRIC  problem  is  solved  by  using  the 
generalized  Lagrange  Multiplier  method  suggested  by  Fox  and  Landi  [Ref.  15:  pp. 
258-261].  Let  O  be  a  Lagrange  Multiplier  associated  with  the  budget  constraint.  The 
Lagrange  function  is  written: 

i     j  i     J 

The  auxiliary  problem  attempts  to  minimize  this  equation.    By  trial 
and  error,  we  try  to  find  the  value  of  O  which  satisfies  a  given  constraint.    Therefore, 
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we  need  to  solve  the  above  equation  for  several  values  of  <t>,  and  choose  that  value  of 
O  for  which  the  required  resources  are  closest  to  the  budget  limit.  Fox  and  Landi  sug- 
gest a  binary  search  procedure.  Their  computational  experience  was  that  at  most  six 
bisections  were  required  to  obtain  budget  allocations  that  were  within  one  half  of  1% 
of  the  original  budget  C. 

The  objective  function  is  separable  in  the  items.    Dropping  the  sub- 
script i  in  the  original  problem  allows  us  to  rewrite  the  equation  as: 


min    >BOi(Sj  I  *JTj)    ~    ®cisu   ~    ®cis, 
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Since  BOj{St\  AjTj)  is  discretely  convex  for  a  given  Si0,  the  optimum 
base  level  is  obtained  by  simply  finding  the  smallest  non-negative  integer  satisfying: 

BO^   +    \\XjTj)    -    BOiSjlljT^QQSio 

(5)  Model  Evaluation.  The  current  Korean  Air  Force  Logistics  system's 
objective  is  maximizing  operational  availability  subject  to  a  budget  constraint.  This  is 
similar  to  the  METRIC  model.  The  objective  of  minimizing  the  expected  number  of 
backorders  is  related  to  the  objectives  of  the  Korean  Air  Force  Logistics  system.  How- 
ever, two  assumptions  are  distinctly  violated.  One  is  batch  size  repair  policy  is  not  used, 
the  other  is  lateral  resupply  is  ignored.  Currently,  these  assumptions  are  not  used  in  the 
Korean  Air  Force  because  of  relatively  long  procurement  leadtimes  and  repair  capability- 
limitations.  It  may  be  possible  for  the  Korean  Air  Force  computer  resources  to  handle 
the  data  requirements  for  the  METRIC  model,  but  it  is  doubtful  for  actual  usage. 
However,  the  METRIC  model  is  obviously  related  to  the  Korean  Air  Force  in  terms  of 
its  basic  concept  and  system  objective. 
b.     Mod-METRIC  Model 

( 1 )  Background.  Mod-METRIC  model  was  developed  by  Muckstadt  to 
deal  with  problems  which  the  METRIC  approach  did  not  consider.  Mod-METRIC 
considers  the  relations  in  parts  hierarchy  and  tries  to  solve  this  multi-indenture  level 
problem.  Mod-METRIC  was  implemented  by  the  U.S.  Air  Force  as  the  method  for 
computing  repairable  stock  levels  for  the  F-15  weapon  system.  [Ref.  16:  p.  472] 

Most  repairables  contain  subassemblies  or  other  components  which 
are  also  repairable.  For  example  an  aircraft  engine,  it  may  have  modules  for  intake, 
combustion  and  exhaust.    If  an  engine  fails,  it  is  replaced  by  a  serviceable  engine  from 
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base  stock.  The  failed  component  is  then  repaired  either  at  the  base  or  depot  level  de- 
pending upon  the  complexity  of  repair.  A  serviceable  module  from  the  base  stock,  if 
available,  will  replace  the  failed  module,  and  the  repaired  engine  is  placed  in  base  engine 
stock. 

(2)  Mathematical  Assumptions.  The  assumptions  stated  in  METRIC  are 
also  applicable  to  Mod-METRIC  except  that  Mod-METRIC  assumes  that  the  demand 
process  is  the  simple  Poisson  process.  [Ref.  16  pp.  475-476] 

(lj  Model  Formulation.  In  METRIC,  the  objective  is  to  minimize  ex- 
pected backorders  for  all  items  subject  to  an  investment  constraint.  The  Mod-METRIC 
objective,  however,  is  to  minimize  the  backorders  for  the  end  item,  subject  to  an  invest- 
ment constraint  on  the  total  dollars  allocated  to  the  end  item  and  its  components  [Ref. 
16.  17:  p.  4~5,  p.  39].  This  difference  is  caused  by  the  hierarchical  maintenance  re- 
lationship between  the  module  and  the  end  item.  As  an  example,  an  engine  backorder 
indicates  that  an  aircraft  is  missing  an  engine  and  is  unavailable  to  perform  its  flying 
mission.  Modules,  on  the  other  hand,  are  used  only  to  repair  engines.  A  backorder  for 
a  module  only  delays  the  repair  of  an  engine.  The  impact  of  module  backorders  and 
engine  backorders  is  clearly  not  the  same. 

Let  7*  denote  the  average  engine  resupply  time  at  a  base.  By  the  na- 
ture of  repairables,  Tt  depends  on  several  factors.  When  it  is  repaired  at  base  level,  T, 
would  be  the  time  it  takes  to  flow  through  the  base  maintenance  system.  If  it  is  repaired 
at  the  depot,  T,  consists  of  the  time  to  place  the  depot  order  for  a  serviceable  part  and 
to  receive  the  part  from  the  depot,  assuming  a  serviceable  asset  is  on  hand  at  the  depot. 
However,  when  there  are  no  serviceable  assets  on  hand  at  the  depot,  an  additional  delay 
is  included  in  the  resupply  time: 

Tt  =  rfr  +  (1  -  rt){At  +  d(S0D)} 

where: 

r,  =  the  probability  an  engine  will  be  repaired  at  base  i 

B,  =  the  average  resupply  time,  given  an  engine  is  repaired  at  base  i 

A,  =  the  average  order  and  ship  time  for  an  engine  at  base  i 

S0  =  the  depot  engine  stock 

D  =  the  average  depot  repair  time 

3S0D  =  the  delay  days  per  demand 
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5S0D  can  be  derived  in  the  following  manner.   The  expected  number 
of  engines  back  ordered  at  the  depot  is: 


BO{S0\/.D)  =    V  (X  -  S0)P(XVD) 


x>s, 


where: 

X  =  Zd     -    r,)A, 

/,  =  the  daily  engine  removal  rate  at  base  i 


In  other  words,  this  expression  is  the  expected  number  of  units  on 
which  delay  is  being  incurred  at  a  random  point  in  time.  Dividing  this  expression  by  the 
expected  number  of  demands  per  day  yields  a  statistic  which  has  the  dimension  of  delay 
days  per  demand: 

SSCD  =  B(S0\aD)Ia 

=  expected  backordcrs'expected  daily  demand 

Further,  B,  can  be  divided  into  two  components.  That  is,  B,  is  equal 
to  the  average  remove  and  replace  time,  given  that  the  necessary  module  is  available, 
plus  the  expected  delay  due  to  the  unavailability  of  the  module  which  is  required  to  re- 
pair the  engine.   Then: 

Bt  =  Rt  +  A,- 

where: 

R,  =  the  average  repair  time  at  base  i  if  modules  are  available 

A,  =  the  average  delay  in  base  engine  repair  due  to  unavailability  of  a  needed  module 

Further,  assume  that  the  engine  should  be  repaired  by  the  failure  of 
module  j.  Then,  A,;  is  the  expected  delay  in  engine  base  repair  time  due  to  a  back  order 
on  module  j  at  base  i.   Thus: 


Au  =   Z  W  ~  S<W  I W  /  h 


x0>so 


where: 
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/.,   =  average  number  of  daily  removals  of  module  j  at  base  i 

S    =  stock  level  of  module  j  at  base  i 

T   =  average  resupply  time  for  module  j  at  base  i 

Then: 

Tij  =  r-fy  +  (1  -  ryXAy  +  djDj) 

where: 

r    =  the  probability  that  a  failure  isolated  to  module  j  will  be  repaired  at  base  level 
B    =  average  base  repair  time  for  module  j  at  base  i 
A    =  average  order  and  ship  time  for  module  j  at  base  i 
D    =  average  depot  repair  time  for  module  j 
8j=    £  (X  -  Sc)P(X\6D)ieD 

X;  S. 

S,   =  the  stock  level  of  module  j  at  the  depot 

=    I-     1    -   r) 
=i 

Consequently,  the  expected  delay  in  engine  repair  at  base  i  due  to 
module  unavailabilitv  is: 


4  =  ( i  /  fti^-iA 


7=1 

Thus,  expression  T,  represents  all  the  system  components  including 
the  depot  engine  and  modules  stock  level,  and  the  base  module  stock  level.  By  using  this 
relationship,  the  engine  and  module  stock  levels  can  be  derived.  Finally,  the  math- 
ematical statement  of  the  problem  is: 


rmn  X  Z  (A^  "  5<)P(A' '  w 


;=1  A">5, 
/  J 


subject  to    V(CESi  +  YjCiJS<l>  +  YjCJS°J  +  CeSo^C 


<=i  j=\  j=\ 


where: 

S  =  stock  level  of  base  i 
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CE  =  unit  cost  of  an  engine  (end  item) 
Cj  =  unit  cost  of  module  j 
C  =  dollar  budget  limit 

(4)  Solution  Technique.  Unfortunately,  the  equation  above  is  not  sepa- 
rable, because  T0  is  a  complex  function  of  the  S,y  Thus,  Muckstadt  recommended  that 
it  should  be  broken  down  into  two  parts;  the  component  subproblem  and  the  end  item 
subproblem  [Ref.  16:  p.  477].  Even  after  the  break-down  of  the  problems,  however,  it 
requires  the  solution  of  many  subproblems  each  of  which  corresponds  to  a  particular 
division  of  the  available  budget  between  components  and  end  items.  The  solution  pro- 
cedure is  as  follows. 

First  allocate  the  budget  C  into  C,  and  C2  for  components  and  end 
items  repectively.  Allocate  C,  among  components  so  as  to  minimize  the  expected  end 
items  repair  delays  summed  over  all  bases  subject  to  the  budget  constraint  C, .  This 
problem  is  mathematically  represented  as: 


mir 

j 

; 

subject  to  2_j(CJSoj  +  Z_jCJS^  ~  Cl 


7=1  /-] 

This  problem  can  be  solved  using  the  METRIC  technique.  Given  the 
result  from  this  step,  compute  the  average  resupply  time  Tt  for  the  components  of  each 
base.  Then,  allocate  the  remaining  budget  C2  so  as  to  minimize  the  expected  end  item 
base  backorders.  The  METRIC  budget  allocation  procedure  may  again  be  used.  The 
steps  provide  a  set  of  proposed  stock  levels  upon  a  given  allocation  of  the  budget  among 
the  end  items  and  components.  These  steps  then  are  repeated  several  times  using  new 
values  for  Q  and  C2  to  establish  the  best  allocation. 

(5)  Model  Evaluation.  Basically,  the  Mod-METRIC  model  has  a  similar 
structure  to  METRIC.  All  assumptions  used  in  Mod-METRIC  model  are  applicable  to 
METRIC  except  for  the  compound  Poisson  demand.  In  other  words,  this  model  is  the 
same  as  METRIC  for  the  Korean  Air  Force  repairable  management  system.  The  only 
disadvantage  is  that  the  Mod-METRIC  model  was  developed  for  new  weapon  systems 
such  as  the  F-15  aircraft.    However,  the  Korean  Air  Force  has  a  lot  of  non-standard 
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weapon  systems.    These  non-standard  weapon  systems  are  causing  problems  such  as 
long  procurement  leadtimes. 

c.     Dyna-METRIC  Model 

(1)  Background.  Dyna-METRIC  was  developed  by  the  RAND  Corpo- 
ration to  provide  an  analytic  method  for  studying  the  transient  behaviour  of 
component-repair  and  inventor}-  systems  under  time-dependent  operational  demands 
and  logistics  decisions  like  those  that  might  be  experienced  in  wartime  [Ref.  14,  18]. 
Note  that  the  past  work  regarding  the  repairable  item  stockage  prior  to  Dyna-METRIC 
only  dealt  with  a  steady  state  inventory  system  with  constant  average  demand  and  ser- 
vice rates.  These  steady  state  assumptions  may  provide  a  good  approximation  during 
peace  time  operations.  But  in  wartime  demands  for  components  may  suddenly  increase 
very  much  relative  to  the  previous  peacetime  operation  and  then  may  decrease  gradually 
or.  in  some  case,  drastically  due  to  attrition  in  the  system.  [Ref.  14:  pp.  1-6] 

A  key  characteristic  of  Dyna-METRIC  is  its  ability  to  deal  with  the 
dynamic  or  transient  demands  placed  on  component  repair  and  inventory  support 
caused  by  time  variables  in  a  scenario  that  includes  sortie  rates,  mission  changes,  phased 
arrival  of  component  repair  resources,  interruptions  of  transportation,  and  the  like,  all 
of  which  would  be  experienced  in  wartime.  It  computes  how  given  resource  levels  and 
process  times  would  contribute  to  the  wartime  capability.  By  exploiting  the  mathemat- 
ical structures  of  its  underlying  equations,  Dyna-METRIC  suggests  the  alternative  cost 
effective  repair  or  stockage  resource  purchases  which  would  achieve  a  target  aircraft 
availability  goal  throughout  the  wartime  scenario. 

Dyna-METRIC  considers  a  three  echelon  inventory'  repair  system. 
Each  base  has  an  in-house  repair  facility  which  may  have  various  test  and  repair  capa- 
bilities. This  base  repair  facility  may  be  supported  by  several  Centralized  Intermediate 
Repair  Facilities  (CIRFs).  Each  operating  base  is  capable  of  doing  only  limited  types 
of  maintenance,  usually  limited  to  simple  removal  and  replacement  operations  at  the 
flight  line.  It  should  be  noticed  that  some  of  the  bases  have  a  CIRF  while  others  do  not 
have  direct  flows  of  parts  to  the  depot.  A  depot  is  represented  as  existing  outside  of  the 
model.  It  is  seen  from  the  model's  point  of  view,  as  an  infinite  source  of  supply  located 
some  order  and  ship  time  away. 

{2j     Mathematical       Assumptions.     The        major        assumptions        of 
Dyna-METRIC  which  distinguish  it  most  from  the  others  are  shown  below. 

1.   Demand  for  items  are  generated  by  a  nonhomogeneous  Poisson  process  with  in- 
tensity function  m(t),  and  mean  value  function  r(t)=   fm{s)  ds.    The  functions, 
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m(t)  and  r(t),  will  be  defined  later.    The  t  denotes  an  arbitrary  time.    Thus,  the  de- 
mand process  is  dynamic. 

2.  The  repair  process  is  independent  of  the  arrival  process  and  has  slack  repair  ca- 
pacity so  that  each  repairable  item  demanding  repair  immediately  receives  services 
with  average  service  time  based  on  the  function  F(s,t).    F(s,t)  will  be  defined  later. 

3.  No  more  than  one  subcomponent  can  fail  or  be  demanded  in  the  repair  of  each 
assembly. 

(3)  Time  Dependent  Pipeline  Equations.  Since  the  major  objective  of  the 
system  is  to  avoid  the  loss  of  aircraft  mission  capability  due  to  a  shortage  of  correctly 
functioning  components  on  the  aircraft,  it  is  necessary  to  compute  the  number  of  com- 
ponents awaiting  repair,  being  repaired,  being  on  the  way  to  and  from  another  echelon 
of  repair,  and  partially  repaired  but  awaiting  spare  parts.  Each  state  is  a  pipeline  seg- 
ment is  which  characterized  by  a  delay  time  that  arriving  components  must  spend  in  the 
pipeline  before  exiting  the  segment.  The  model  expands  each  component's  expected 
pipeline  size  into  a  complete  probability  distribution  for  the  number  of  components 
currently  undergoing  repair  and  on  order,  so  the  probability  distribution  for  all  compo- 
nents can  be  combined  to  estimate  aircraft  availability  and  sorties.  [Ref.  14:  pp.  7-23] 

Under  the  assumptions  that  the  probability  distribution  of  repair 
time  is  independent  of  the  failure  process,  the  average  number  of  components  in  the  re- 
pair pipeline  will  be: 

ka  =  dT 

where: 

d  =  average  daily  failure  rate 
7"  =  average  repair  time 

With  the  further  assumption  that  demand  has  a  Poisson  probability  distribution,  the 
probability  that  there  are  K  components  in  the  pipeline  at  any  point  in  time  would  be: 

P{  K  in  pipeline  )  =  ).fse~ }"  /  K ! 
However,  in  Dyna-METRIC  the  demand  is  a  function  of  time  so  that: 
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d(t)  =  (failure  per  flying  hour) 
x( flying  hour  sorties  at  time  t) 
x(number  of  sorties  day  per  aircraft  at  time  t) 
x(quantity  of  the  component  on  the  aircraft) 
x(percentage  of  aircraft  with  the  component) 

Also,  in  place  of  a  constant  average  repair  time,  T,  the  dynamic  model  uses  the  proba- 
bility that  a  repair  started  at  time  s  is  not  completed  at  time  t.   That  is: 
F(t.s)  =  Probability(component  entering  at  s  is  still  in  repair  at  t) 
=  Probability( repair  time  >  t-s  when  started  at  s) 

The  average  number  of  components  in  the  pipeline  is  derived  by  combining  those  two 
functions.  Consider  only  those  components  that  arrived  in  an  interval  of  time.  As.  cen- 
tered at  time  s.   Then: 

A/.(r,s)  =  d(s)  F{i,s)As 

where: 

A/.(/.5)  =  expected  number  of  components  in  the  repair  pipeline 
dis)  =  daily  failure  rate  at  time  s 

F{i.s)  =  Probability  that  a  component  is  not  out  of  repair  by  time  t 
As  =  interval  of  time  centered  at  s 

If  we  assume  that  the  number  of  failures  arriving  in  the  interval  As  is  independent  of  the 
number  of  failures  arriving  in  similar  intervals  centered  at  other  time  other  than  s  and 
Fit)  is  independent  of  the  probability  distribution  generating  the  demand  rate,  then: 

A(r)  =  £y(5)  F(t,s)  As 
Further  assume  that  As  is  very  small,  so  that: 


/(/)  =  ^'d(s)F(t,s)ds 
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With  the  additional  assumption  that  the  component  failure  probabil- 
ity distribution  is  Poisson,  ).{t)  is  the  mean  of  a  time  varying  (nonhomogeneous)  Poisson 
process.   That  is,  the  probability  of  K  components  in  repair  at  time  t  is: 

P{K)  =  X{tfe~mlK\ 

where: 


;.(/)  =  j'd{s)F{t,s)ds 


(4)     Time  Dependent  Component  Performance  Measure.      The  component 
measures  typically  computed  by  the  Dyna-METRIC  model  are  [Ref  14:  pp.  24-29]: 
R(t)  =  ready  rate  at  time  t 
FR(t)=  fill  rate  at  time  t 
EB(t)=  expected  back  orders 
VBD(t)  =  variance  of  the  backorders 
DT(t)=  average  cumulative  demands  by  time  t 
S(t)=  supply  level  at  time  t 


The  readv  rate  is  eiven  bv: 


5(r) 

R{t)  =  y  p{{K\).{t))} 


A=0 


Since  the  definition  of  the  fill  rate  is  the  probability  that  a  component 
will  be  available  when  a  demand  is  placed,  it  is  therefore  the  probability  that  demands 
have  left  at  least  one  component  available,  that  is,  the  sum  of  the  probabilities  of  de- 
mands less  than  the  stock  level.   Expected  backorders  are  given  by: 


EB{t)  =    £{*-  s(t)}P{K  /;.(/)} 


K>S(t) 


=  ;.(/)  -  s(t)  +  £(0{*(0  -  K}P{Kj).{t)} 


A'=0 
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For  K  greater  than  s(t),  there  will  be  backorders  of  {K  -  s{t)}.  The 
probability  of  any  demand  level,  that  is  K,  is  P{K /  /(/)},  and  the  expected  value  of  the 
backorders  is  merely  the  product  of  the  various  values  the  backorders  can  take  on  times 
the  probability  of  a  demand  at  that  given  value.  The  variance  in  backorders  is  given  by: 

VB(t)  =    ]T  iK  ~  s(t)?P{Kim  ~  imt)}2 

K>S(r) 

(5)  Time  Dependent  Optimal  Determination  of  Spare  Parts.  The  fact 
that  pipelines  have  time-dependent  probability  distributions  means  that  the  optimal  mix 
of  spare  components  at  one  point  in  time  may  not  be  the  optimal  mix  at  another.  Thus, 
the  approach  to  take  is  to  compute,  for  each  time  period  of  interest,  the  marginal  in- 
crease in  spare  parts  to  achieve  a  given  capability  over  those  already  input  or  determined 
for  a  previous  time.  In  determining  the  supply  level,  the  model  attempts  to  provide 
enough  spare  parts  to  give  the  desired  confidence  at  the  lowest  cost  at  each  point  in  time 
of  interest.   Thus,  the  objective  function  is  to  minimize  the  total  cost  of  spare  parts. 

Let: 

S,  =  the  spare  parts  level  for  component  i 

C,  =  the  unit  cost  of  component  i 

a  =  the  desired  confidence  level 

A',  =  the  non-mission  capable  rate  not  to  be  exceeded 

P{KnS)    =    the   probability   that   the   non-mission   capable   rate   is   less   than   K„ 
given  a  stock  level  S 

I  =  types  of  repairables  required  on  each  aircraft 
Then,  the  problem  to  solve  is: 

/ 


minimize  /  CjSj 


subject  to  P{Kn,S)  >  a,      Sj  >  Sl0 

where: 

S,c  =  the  input  stock  level  or  previous  time  optimization  stock  level  for  component  i 
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Assuming  complete  cannibalization,  P(Kn,S)  equals: 

p(Kn,s)  =  YlnQ^) 

where: 

P'(Q,Kn)  =    z  PM 

K=0 

P,(K)  =  the  probability  of  exactly  K  failures  of  component  i 

The  necessary  condition  for  the  performance  constraint  to  be  met  is  to  have, 
P'(QiK)  ^  <*,  for  each  i.  Then  marginal  analysis  is  used  to  determine  the  best  mix  of 
additional  components  to  achieve  the  desired  goal.  This  process  proceeds  by  investing 
in  one  additional  component  at  a  time  which  is  selected  by  finding  the  component  that 
gives  the  largest  increase  in  the  logarithms  of  the  confidence  level  at  the  lowest  cost. 
That  is,  we  determine: 

am  p  (Kn,S)  i  q 

where: 

MnP(Kn,S)  =  ln{P(Kn,S>)IP{Kn,S)) 

O'     =     (j[,  J),  —  >-J,~n j 

The  component  for  which  supply  is  increased  one  unit  is  the  one  whose  index  solves; 
max  A,  In  P(K„,S)  /  C,.  This  process  continues  until  the  given  confidence  level  is  achieved. 
At  this  point,  the  resulting  value  of  S  is  the  efficient  solution  of  the  base  stockage 
problem. 

(6)  Model  Evaluation.  The  Dyna-METRIC  model  is  relatively  hard  to 
adapt  for  the  Korean  Air  Force.  Especially  since  this  model  assumes  that  each  repair- 
able item  needing  repair  immediately  receives  service.  This  does  not  occur  in  the  Korean 
Air  Force.  Currently,  the  Korean  Air  Force  applies  a  batch  repair  policy  for  repairables, 
except  in  case  of  NORS  condition  aircraft.  This  is  a  serious  violation  of  the  assumption 
of  Dyna-METRIC.  So,  the  Dyna-METRIC  model  is  difficult  to  use  for  the  Korean  Air 
Force  repairable  management  system  without  using  a  continuous  repair  policy. 
d.     Mathematical  Models  Used  In  The  Navy 

( 1 )     Background.      The  Uniform  Inventory  Control  Program  (LTCP)  was 
developed  in  August  1965  to  provide  a  standard  system  to  be  used  at  all  L'.S  Naval 
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Supply  Systems  Command  (NAVSUP)  ICPs.  The  Fleet  Material  Support  Office 
(FMSO)  under  the  direction  of  NAVSUP  is  responsible  for  the  system  design,  ADP 
analysis,  programming  and  documentation  of  LTCP.  The  development  of  the  L'ICP 
formulas  for  inventor}'  levels  follows  the  approach  used  by  Hadley  and  Whitin.  [Ref. 
13.  19] 

(2)  Mathematical  Assumptions 

1.  It  is  a  continuous  review  system.  Wholesale  inventory  levels  requirements  and  as- 
sets are  known  by  the  Inventory  Control  Point  (ICP)  at  all  times.  [Ref.  13,  19:  pp. 
162-165,  Chapter  3  Appedix  A] 

2.  It  is  a  steady  state  environment.  The  key  characteristics  of  the  items  managed  by 
the  ICP  are  assumed  to  be  constant  over  the  forecast  period.  These  characteristics 
include  the  forecasted  average  values  and  variances  of  the  random  variables  of  the 
rate  of  customer  demand,  procurement  leadtime.  production  leadtime,  depot  repair 
times,  depot  repair  survival  rate  and  the  rate  of  carcass  returns. 

3.  An  order  for  procurement  or  repair  is  placed  when  the  assets  reach  the  reorder  level 
or  the  repair  level. 

4.  Customer  demands  and  carcass  returns  do  not  occur  in  more  than  one  unit  per 
transaction. 

5.  The  unit  procurement  cost  or  repair  cost  of  an  item  is  independent  of  the  magni- 
tude of  order  quantity  or  repair  quantity. 

6.  The  cost  of  a  backorder  and  the  time-weighted  cost  of  a  backorder  can  be  accu- 
rately quantified. 

7.  The  reorder  level  and  repair  level  are  always  non-negative. 

8.  The  cost  to  hold  one  unit  of  stock  in  the  inventor}'  is  proportional  to  the  unit  cost 
of  the  item. 

9.  No  interaction  exists  among  families  of  items  or  individual  nonfamily  items  or 
both.  Each  family  or  nonfamily  item's  inventory  levels  requirements  are  calculated 
independently  of  those  of  families  or  nonfamily  items. 

10.  The  optimal  inventory  levels  are  determined  by  minimizing  an  average  annual 
variable  costs  equation  composed  of  order  cost  plus  holding  costs  plus  shortage 
costs. 

(3)  L'ICP  Depot  Level  Repairables  (DLRs)  Procurement  Model.  The 
repairables  procurement  model  starts  with  a  total  variable  cost  (TVC)  equation  which 
is  to  be  minimized.  A  notable  difference  from  the  model  for  consumables  is  the  inclusion 
of  receipt  of  RFI  assets  from  a  repair  process  in  the  DLR  model. 

(a)     Ordering  Cost 


Let: 

A  =  administrative  cost  to  order 
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D  =  forecasted  quarterly  recurring  demand 
B  =  forecasted  quarterly  regenerations 
Q  =  economic  order  quantity 

Then,  A(D  —  B)jQ  is  the  expected  number  of  procurements  per  year.    Thus,  the  annual 
ordering  cost  is  given  by  A(D  —  B)AjO. 
(b)     Holding  Cost 

Let: 

I  =  percent  cost  per  dollar  of  inventory  held  annually 
C  =  item  cost  (replacement) 
L  =  procurement  leadtime 
T  =  repair  turnaround  time 
R  =  repair  level 

F(x)  =  the  function  of  probability  distribution  of  leadtime 
demand 

The  mean  demand  during  a  procurement  leadtime  is  given  by: 

DT  +  {D  -  B)  x  (L  -  T) 
Thus,  holding  cost  is  given  by: 

O  f 

IC\_~y   +  R  -  {{D  -  B)  x  L  +  B  x   7}  +  {x  -  R)F{x)dx'] 

x>R 

let: 

).  =  shortage  cost  per  requisition  short 
E  =  military  essentiality  weight 
F  =  quarterly  requisition  rate 

Then: 


(;.£){4(Z)  -  B)  I  Q)F 
d\     {x  -  R)  F(x)  dx 


'x>R 
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It  equals  the  cost  of  the  expected  number  of  backorders  in  a  year.    The  TVC  equation 
is  symbolized  by: 

TVC  =  (D  -  B)AIQ 

+  ICIQI2  +  R  -  {(D  -  B)xL  +  BxT)  +  f     (x  -  R)F(x)  dx] 

x>R 

+  (>.£)[4(D  -  B)IQ]FIDf     (x  -  R)F(x)  dx 

*>R 

Then,  by  setting  dTVC  /  dQ  =  0: 

Q  =     /[8(Z)  -  B)(A  +  ).EFj  D)\     {x  -  R)F{x)  dx  I IC] 

*  x>R 

Also,  by  setting  dTVC  ,  dR  =  0,  we  find: 

OICD 


i 


(x  -  R)F{x)  dx  = 


{QICD    +    4aEF(D  -  B)} 

x>R 

Because  the  expression  f  (x  —  R)F(x)  dx  is  the  cumulative 
distribution  for  leadtime  demand,  this  is  the  quantity  defined  as  RISK.  Note  that  since 
Q  and  R  are  related,  the  reorder  quantity  Q  cannot  be  solved  independently.  Thus, 
LTCP  approximates  Q  by  using  a  variation  of  the  economic  order  quantity  formula: 


8(Z>  -  B)A 


IC 


Then,  OICD  /  {OICD    +    A).EF{D  -  B)}  can  be  computed  for  RISK  determination. 

f4j  UICP  Depot  Level  Repairables  Repair  Model.  The  repair  model  also 
starts  with  a  total  variable  cost  equation  viewed  largely  independently  of  the  procure- 
ment problem.  Note  that  its  time  horizon  is  the  depot  level  repair  turn  around  time. 
The  total  variable  cost  equation  for  the  repair  model  is  the  sum  of  the  order  cost,  holding 
cost  and  backorder  cost. 
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Let: 

Q2  =  economic  repair  quantity 

A2  =  repair  administrative  order  cost 

C2  =  repair  price 

R2  =  repair  level 

Fz(x)  =  probability  distribution  of  demand  during  during  repair  turnaround  time 

>.2  =  repair  shortage  cost 

Then: 

JVC  =  A  mm(D,A)A2  /  Q2 

+  IC2{-Q-   +  R2-  (DXT)  +  f      (x  -  R2)F2(x)dx) 

£  *>R2 


x>  R-) 


+  /2£{4  min{D,B)  /  Q2)F  /  Df    _  (x  ■  ■  R2)F2{x)dx 
By  setting  8TVC  /  802  =  0: 


0,  =      Umin{D,B)A2  +  >.2EF j  d\       {x  -  R1)F1(x)dx /  IC2 


^  n2 


Also,  by  setting   dTVC  /  6R2  =  0: 

j      F2{x)dx  =  Q2IC2D  I  Q2IC2D  +  A).2EF  min{D,B) 

TOO 

Again,  Q2  and  R2  are  related.   The  LTCP  model  approximates  Q2  as  follows: 


Q2    =    ^8  min  (D,B)A2  /  IC2 
Then,  Q2IC2D  /  QJC2D  +  A).2EFB  can  be  computed  by  using  the  above  results. 
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f5;  Integrated  Repairables  Model.  As  stated  earlier,  the  requirements 
computed  by  the  procurement  and  repair  models  are  accomplished  independently  of 
each  other.  As  a  result,  this  leads  to  a  carcass  constrained  situation.  That  is,  the  com- 
puted procurement  inventor}-  level  for  an  item  does  not  provide  sufficient  carcasses  to 
allow  repairs  at  the  computed  repair  inventory  level.  To  solve  this  problem,  Naval 
Supply  System  Command  (NAVSUP)  has  made  some  changes.  Under  the  model  inte- 
gration, there  is  only  one  RISK  formula: 


RISK  =    IC3D  I IC3D  +  aFE 


where: 


C3    =    (BID)(C2)  +  (1  -    -£-)(Q 

Also,  rather  than  using  a  procurement  leadtime  or  a  depot  level  re- 
pair turnaround  time,  it  is  uses  an  average  acquisition  time  as  the  horizon  for  computing 
the  safety  level.   The  average  acquisition  time  (L2)  is  defined  by: 

L2    =    (1  -  B\D)L    +    {BjD)T 

It  is  clear  that  L2  is  the  weighted  average  of  the  procurement  leadtime 
and  the  repair  turnaround  time  because  D-B  represents  the  quantity  to  be  procured  and 
B  represents  the  quantity  from  the  regeneration. 

(6j  Model  Evaluation.  For  the  LTCP  model,  three  major  problems  exist. 
The  first  problem  is  that  LTCP  models  is  a  continuous  review  inventor}'  system.  The 
Korean  Air  Force  repairable  management  system  applies  a  periodic  review.  The  second 
problem  is  that  LTCP  assumes  that  no  interaction  exists  among  failures  of  items. 
Finally,  the  optimal  inventory  level  for  LTCP  model  is  determined  by  minimizing  annual 
variable  cost  (order  costs  plus  handling  costs  plus  shortage  costs).  However,  the  Korean 
Air  Force  does  not  have  t,he  same  objective. 

e.     Availability  Centered  Inventory  Model 

(1)  Background.  ACI.M  was  developed  by  CACI  Inc.  under  the  spon- 
sorship of  the  Ship  Support  Improvement  Project,  Naval  Sea  System  Command, 
PMS-306,  and  aproved  in  March  1981  by  the  Chief  of  Naval  Operations  for  use  in  de- 
termining consumer  level  stockage  quantities  for  selected  equipment  [Ref.  20:  pp.  1-2]. 
However,  this  model  was  initially  used  as  a  part  of  a  larger  model  (Logistics  Support 
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Economic  Evaluation  Model)  which  provided  necessary  inputs  and  enabled  comparisons 
with  other  Navy  stockage  policies.  After  this,  successive  refinements  for  the  simplifi- 
cation of  the  solution  procedure  and  associated  computer  programs  have  been  made  for 
several  years. 

This  model  was  originally  designed  to  calculate  inventory  levels  for 
all  items  in  the  parts  breakdown  of  an  equipment  and  at  all  stockage  facilities  in  a 
multi-echelon  support  system.  Thus,  ACIM  is  capable  of  computing  levels  for  all  ships, 
intermediate  maintenance  activities,  and  depots  that  use  or  support  the  equipment. 
However,  this  model  is  mainly  used  in  the  provisioning  process  to  compute  shipboard 
allowances  to  achieve  specified  weapon  system  readiness  levels  which  have  not  been 
achieved  with  the  standard  protection  level  models. 

(2)  Axailablity  Measure.  ACIM  recognizes  that  the  purpose  of  a  supply 
system  is  to  provide  sufficient  support  so  that  a  weapon  system  is  operational  when  it 
is  needed.  The  terminology  used  to  described  this  goal  is  operational  availability  (A0). 
ACIM  defines  A„  by  the  following  formula  [Ref.  21:  p.  5]: 


A0  =  up  time  /  (up  time  +  down  time) 

=  MTBF  /  (MTBF  +  MTTR  +  MSRT) 


A0  may  also  be  interpreted  as  the  probability  that  the  equipment  is  in  an  operable  con- 
dition at  a  random  point  in  time.  Among  the  three  factors  of  Ae,  the  MTBF  and  MTTR 
are  system  parameters  outside  the  control  of  the  supply  system  and  are  viewed  as  con- 
straints. MSRT  is  the  only  term  which  depends  on  stockage  postures,  it  is  therefore  the 
one  that  the  ACIM  model  focuses  on  to  achieve  a  given  value  for  A0. 
(3)     Mathematical  Assumptions 

1.  All  parts  are  organized  in  terms  of  equipment  with  a  top-down  breakdown  that  can 
be  represented  as  an  arborescent  network.  Any  part  may  be  totally  consumable, 
totally  repairable,  or  any  mix  thereof.  [Ref.  21:  pp.  10-11] 

2.  Stockage  and  maintenance  facilities  are  organized  in  a  hierarchical  structure  ac- 
cording to  supply/  maintenance  flows  which  can  be  represented  as  an  arborescent 
network  as  illustrated  below.  Each  facility  has  a  colocated  maintenance  and  supply 
capability.  The  facility  at  the  top  of  the  structure  is  assumed  to  have  an  infinite 
supply  of  all  items. 

3.  External  demands  upon  supply  are  stationary  and  compound-Poisson  distributed. 

4.  All  stockage  locations  use  a  continuous  review,  (S-1,S)  ordering  policy. 
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5.  Mean  Time  to  Repair  (MTTR)  is  defined  to  include  all  equipment  downtimes  that 
are  not  supply  related. 

6.  There  is  no  lateral  resupply  (transshipment)  among  bases. 

7.  Items  repaired  at  any  location  are  assumed  to  be  returned  to  collocated  stocks  for 
issue. 

S.  If.  for  a  given  facility  in  the  network,  the  stocks  are  physically  distributed  in  several 
places,  it  is  assumed  that  the  resupply  time  for  direct  customers  (next  lower  ech- 
elon) is  independent  of  the  location. 

9.  The  order  policy  assumption  precludes  consideration  of  economies  of  scale  for  re- 
supply.  That  is,  all  ordering  is  on  a  one  for  one  basis. 

(4)  Model  formulation.  The  goal  of  AC1M  is  to  maximize  the  opera- 
tional availability  (Ae)  of  a  weapon  system  subject  to  a  given  inventory  budget.  With  the 
definition  given  above  of  A,  and  with  the  view  that  MSRT  is  the  only  term  which  is  af- 
fected by  the  stockage  decision,  the  developers  of  ACIM  argue  that  the  allocation  which 
maximizes  Ac  is  equivalent  to  the  allocation  which  minimizes  MSRT.  ACIM  therefore 
actually  attempts  to  minimize  MSRT  subject  to  given  constraints. 

Let  i  be  an  arbitrary  item  in  equipment  e  (which  may  be  e  itself).  Let 
u=0  represent  an  arbitrary  facility  in  the  support  system  and  u=  1,2,3,....,  represent  fa- 
cilities at  the  next  lower  level,  i.e.  those  facilities  that  submit  items  for  repair  directly  to 
or  obtain  resupply  from  facility  o.  Then,  the  objective  can  be  explicitly  stated  as  follows: 


Find  values  for  Skv  which  minimize  D,u  for  all  user  locations  u: 


subject  to  yCkSkv  =  B 


fc.v 

where: 

Ck  =  unit  cost  of  item  k 

B  =  given  budget  for  spares  procurement 

D,u  =  expected  delay  per  demand  upon  inventory  for  for  item  i  at  location  u 

Skt  =  the  level  of  inventory  for  item  k  at  location  v 

DIU  is  one  of  the  components  of  Mlu  which  is  defined  as  the  mean  time  to  return  a  failed 
unit  of  item  i  at  location  u  to  a  serviceable  condition.   M,u  is  given  by: 
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where: 

D,u  =  expected  delay  per  demand  upon  inventory  for  for  item  i  at  location  u 
Tlu  =  mean  time  to  repair  item  i  at  user  location  u 

Diu  in  turn  is  given  by: 


Dlu  =  (1/4J^(A'-  Slu)P{X:?.iuTiu) 


where: 

Sw  =  stock  level  of  item  i  at  location  u 

/.,u  =  expected  number  of  demands  upon  inventory  for  item  i  at  location  u 
P(X:  /.,uTlu)  =  probability  of  X  units  of  stock  reduction  for  item  i  at  location  u 
Ttu  =  mean  resupply  time  for  item  i  at  location  u 

Tiu  is  given  by: 

Tiu  =  (Llu  +  £,*)+(!-  rlu)(Rlu  +  R-u) 

where: 

r,u  =  probability  that  a  demand  for  item  i  upon  inventory  at  location  u 

Llu  =  average  resupply  leadtime  assuming  stock  is  is  available 

L'u  =  additional  resupply  leadtime  due  to  expected  shortage  at  the  resupply  source 

Rlu  =  average  shop  repair  cycle  assuming  availability  of  spares  for  items  within  i 

R'u  =  additional  shop  repair  cycle  due  to  expected  shortage  of  spares  items  within  i 

Since  ACIM  assumes  that  MTBF  and  MTTR  are  independent  of  the 
stockage  policy,  minimizing  D,u  is  assumed  to  maximize  Aeu  (operational  availability  of 
a  equipment  at  location  u). 

(5)  Solution  Technique.  Assume  that  initial  values  for  SIU  are  given  for 
all  items  and  locations.  These  may  all  be  zero  or  some  minimum  value  given  by  policy 
or  current  assets.   Then  compute  the  MSRTs  for  all  items. 

ACIM  calculates  what  the  new  MSRT  would  be  if  one  additional 
unit  of  stock  were  placed  against  that  particular  item.  Subtracting  the  new  MSRT  from 
the  old  MSRT  and  then  dividing  by  the  unit  cost  of  the  item  provides  the  model  with  a 
value  which  is  multiplied  by  the  item's  demand  to  determine  a  selection-rank.  After  the 
selection  number  is  calculated  for  each  item,  one  unit  of  stock  is  added  to  the  candidate 
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with  the  highest  selection  rank  number.    The  operational  availability  (A„)  is  calculated 
and  if  the  target  Ag  has  not  been  achieved,  the  model  continues  the  same  procedure. 

f6)  Model  Evaluation.  As  with  the  UICP  model,  three  major  violations 
of  the  basic  assumptions  of  the  ACIM  model  exist.  The  first  is  that  ACIM  assumes  that 
you  have  a  continuous  review  system.  The  second  is  that  ACIM  assumes  there  is  no 
lateral  resupply.  Finally,  ACIM  assumes  that  the  stocks  are  physically  distributed  at 
several  locations.  Obviously,  the  Korean  Air  Force  supply  system  does  not  require  se- 
veral physical  locations,  even  though  the  system  objective  is  to  maximize  operational 
availability. 

D.     SUMMARY 

Throughout  Section  B  and  C  of  this  chapter,  we  have  reviewed  the  various  inven- 
tory models  and  theories  for  repairable  items  which  are  currently  in  use  in  the  military 
services,  or  in  the  process  of  being  evaluated.  All  of  the  assumptions  and  the  math- 
ematical formulations  of  the  models  were  stated  and  discussed. 
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IV.     MODEL  PROPOSAL  FOR  THE  KOREAN  AIR  FORCE 

A.  INTRODUCTION 

In  Chapter  II,  we  summarized  the  current  Korean  Air  Force  repairables  manage- 
ment system  and  repair  process.  The  general  structure  of  repairables  management  in  the 
Korean  Air  Force  is  similar  to  the  structure  assumed  by  several  of  the  models  described 
in  Chapter  III. 

In  this  Chapter,  we  describe  a  proposed  inventory  control  process  for  the  Korean 
Air  Force.  This  process  is  based  partly  on  queueing  models,  and  partly  the  Wilson- 
Harris  economic  order  quantity  model. 

B.  THEORETICAL  BASE  FOR  THE  MODEL  PROPOSAL 

1.  Notation 

c  =  number  of  identical  repair  channels 

ix,  =  average  repair  rate  per  repair  channel 

a,  =  average  failure  rate  of  an  individual  unit  of  a  repairable  item 

X,  =  random  variable  describing  the  steady 
state  number  of  NRFI  units  of  item  i 

P„  =  steady  state  probability  that  there  are  n  NRFI  units 

2.  The  M/M/l/K/K  Queueing  System 

This  model  is  a  limited  source  queueing  model  in  which  there  are  only  K  cus- 
tomers. It  is  variously  called  the  machine  repair  model,  the  machine  interference  model, 
or  the  cyclic  queueing  model.  It  is  one  of  the  most  useful  of  all  queueing  theory  models. 
One  way  to  view  this  model  is  shown  in  Figure  6.  [Ref.  22:   pp.  186-190] 

The  population  of  potential  customers  for  this  queueing  system  consists  of  K 
identical  devices,  each  of  which  has  an  operating  time  of  O  time  units  between  break- 
downs, O  having  an  exponential  distribution  with  average  value  1  /  a.  The  one 
repairman  repairs  the  machines  at  an  exponential  rate  with  an  average  repair  time  of 
1  /  ix  time  units.  The  operating  machines  are  outside  the  queueing  system  (outlined  by 
the  dashed  lines)  and  enter  the  system  only  when  they  break  down.  The  queueing  system 
always  reaches  a  steady  state  because  there  can  be  no  more  than  K  customers  in  the 
system  (one  machine  being  repaired  and  K-l  waiting  for  repairs). 
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Figure  6.      Queueing  System  with  One  Repairman  (M/M/l/K/K) 

When  n  of  the  machines  arc  down  (not  operating),  then  K-n  of  them  arc  oper- 
ating and  the  time  until  the  next  machine  breaks  down  is  the  minimum  of  K-n  identical 
exponential  distributions  with  parameter  (K-n)  a. 
3.     The  M/M/c/K/K/  Queueing  System 

This  queueing  system  is  similar  to  the  machine  repair  model  considered  in  the 
section  above  except  that  we  have  c  rather  than  one  repairman,  where  c  <  K,  as  shown 
in  Figure  7.  [Rcf.  22:  pp.  190-192] 

C.     DEVELOPING  A  MODEL  FOR  THE  KOREAN  AIR  FORCE 
1.     Application  Conditions 

1.  Supporting  D  items,  i.e..  we  want  D  items  to  be  available  at  all  times.  In  this  case, 
D  represents  the  number  of  aircraft  which  the  Korean  Air  Force  wants  to  have 
available  at  all  times. 

2.  P(Succcssful  Repair)  is  close  to  1.  When  P(successful  Repair)  =  1,  we  have  a  pure 
machine  repair  queueing  system.  When  P(SucccssfuI  Repair)  =  0.  we  have  a  con- 
sumable item  whose  management  is  best  accomplished  using  standard  methodol- 
ogies from  inventory  theory.  The  sample  data  collected  from  the  Korean  Air  Force 
for  this  study  indicated  that  P(Succcssful  Repair)  =  0.9S  or  higher  for  all  repairable 
items  examined.  For  this  reason,  the  proposed  model  is  based  partly  on  the  ma- 
chine repair  queueing  model.   The  data  base  used  for  the  study  was  obtained  from 
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Figure  7.      Queueing  System  with  c  Repairman  (M/M/c/K/K) 

the  KAFLC  during  the  19S4-1988.    Table  3  shows  the  mean  probability  of  suc- 
cessful repair  based  on  the  sample  data. 6 

Table  3.      P(SUCCESSFUL  REPAIR)  COMPUTED  FROM  SAMPLE  DATA 


Year 
AC 

84 

8.5 

86 

87 

88 

Average 

F-16 

NA 

NA 

1.000 

l.ono 

1.000 

1.000 

\:-5 

0.98  S 

0.980 

0.987 

0.986 

0.984 

0.985 

I  :-4 

0.986 

0.987 

0.986 

0.976 

0.987 

0.9S4 

3.  Conditions  close  to  the  M/MT.K/K  queueing  system  hold.  Because,  there  are  a 
number  of  repairmen  in  repair  depot,  it  might  be  reasonable  at  first  to  assume  that 
the  M/M/c/K/K  model  applies.  However,  most  repairs  require  specialized  skills 
held  by  few  people  for  a  given  repairable  part.  Further,  the  repair  depot  is  stalled 
and  scheduled  so  that  there  is  little  overall  excess  capacity.  As  a  result,  there  is 
little  change  in  the  repair  process  output  rate  when  the  second,  third,  fourth,  etc. 
unit  backs  up  in  the  repair  queue. 

4.  Steadv  state  has  been  achieved. 


6   The  sample  data  list  appears  in  Appendix  A  and  B 
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2.     Model  Formulation 

This  proposed  model  consists  of  3  steps.  The  first  step  computes  a  minimum 
population  size  A',  for  the  ith  type  of  repairable  item.  Population  refers  in  this  case  to 
the  installed  material  plus  the  stored  RFI  material  plus  the  NRFI  material  awaiting  re- 
pair.   In  other  words,  all  material,  both  good  and  broken. 

Once  the  minimum  population  size  K,  is  chosen,  an  attrition  buy  safety  level 
S,  is  chosen.  This  safety  level  is  used  to  insure  that  the  actual  population  size  usually 
stays  above  K,  during  the  period  when  the  system  is  awaiting  the  arrival  of  an  attrition 
buy. 

The  last  step  is  to  compute  the  size  of  an  attrition  buy  itself.    This  Q]  is  com- 
puted using  the  Wilson-Harris  economic  order  quantity  model. 

In  step  I,  described  on  the  next  page,  the  M  \1  1/K/K  queueing  system  will 
be  used  to  obtain  probability  expressions  of  the  form  P{S\>  A'  —  D).  These  expressions 
are  obtained  as  follows: 

k, 
P(A-;>A^-Z)()  =         Yj       W^) 

n=  K,  -  D,  +\ 

where: 

.V  =  Discrete  random  variable  representing  the  number 
of  units  of  item  i  which  are  broken 

A  =  Decision  variable  representing  planned  minimum 

number  of  units  of  the  item  to  be  owned,  including  those 
installed  in  the  repairable  equipment 

D  =  Number  of  units  of  item  i  which  the  Korean  Air 
Force  wants  to  have  available  at  all  times 

The  right  hand  side  of  the  last  expression  is  equal  to: 

k,  -  D, 

1  -    J]  pw  =  ») 

using  the  M  M  IKK  queue,  Allen  [Ref.  22:  p.  192]  gives  the  following  formula  for 
P(A>/7): 
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where: 


fc=0 

In  step  2,  described  in  the  following  section,  the  expression  P(X:  >  S,)  is  also 
used.    It  is  obtained  using  the  Poisson  distribution: 

{DaiLife-{D°iL') 
P(Xt  =  K)  =   — 


A! 


where: 


X,  =  discrete  random  variable  representing  attrition 

demand  for  item  i  during  a  procurement  leadtime 

Dai  =  forecasted  attrition  demand  for  item  i 

A,  =  forecasted  procurement  leadtime  for  item  i 

We  will  want  S,  such  that  P(X,  >  S,)  is  less  than  a  risk  value,  b2„  which  is  pre-specified 
by  the  user,  where  P{Xt  >  S,),  is  given  by: 

PiXt>S}  =     Y,   pW  =  k"> 

k=S,+\ 


-    I 


(4,A)V(D-« 


k\ 

k=S.  +1 


i  _  v  (DaiLfe-^ 


a.     Step  I 

There  are  several  different  objective  functions  which  might  be  used  to  help 
choose  the  size  of  the  population  fot  the  ith  item.  The  following  four  subsections  de- 
scribe four  different  math  programs  for  choosing  K,. 

( 1)     Alternative  I 

Choose  A'  by  the  following: 


58 


min     Kt 
subject  to  P(Arj  >  {Kt  -  £>,))  <  bu 


where: 


D  =  Number  of  units  of  item  i  which  the  Korean  Air 
Force  wants  to  have  available  at  all  times 

b1:  =  Parameter  chosen  by  the  decision  maker 

This  first  alternative  formulation  contains  no  explicit  budget  constraint.  However, 
minimizing  A'  for  all  items  is  equivalent  to  minimizing  the  cost  of  that  portion  of  the 
population  which  is  owned  to  meet  the  availability  goals.  Unfortunately,  this  formu- 
lation doesn't  identify  solutions  which  exceed  the  available  budget.  The  next  math 
programming  formulation  accomplishes  this. 
r2j     Alternative  II 


m 

min 


2> 

subject  to  P(X;  >  Kt  —  D^)<bu,        i=  1, ,m 

m  m 

and    /A^i  +  ^fl/[max(0,  Kx  -  Koi)\  <  B 

where: 

h,  =  annual  holding  cost  per  unit  held 

B  =  total  budget  for  one  year 

m  =  number  of  repairable  items  being  considered 

K0i  =  previous  planned  minimum  number  of  units  owned 

a,  -  unit  cost  of  item  i 

r3j     Alternative  III 
Maximize  customer  service  subject  to  investment  constraint: 
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m 


z  =  ^]p(a>(a;-  a)) 


subject  to  /  hjKj  +  )>  «,[max(0,  Kt  -  Kol)]  <  B 

(4)     Alternative  IV 
Maximize  customer  service  subject  to  investment  constraint  (another  form): 

min    Z  =  max  {/>(#,>(£)  -  Z)())} 

m  m 

subject  to  /  hjKj  +  /^ a,-[max(0,  Kt  -  Kol)]  <  B 

The  choice  between  these  alternate  formulations  depends  on  the  way  the  Korean  Air 
Force  feels  most  comfortable  with  resolving  the  conflict  between  budget  limitations  and 
customer  service.  For  the  numerical  computation  of  this  study,  we  chose  alternative  II. 
This  choice  balanced  the  need  for  realism  against  the  time  constraints  imposed  on  this 
research. 

b.     Step  2 

Attrition  buy  reorder  point,  R,  =  K,  +  St,  where  S,  is  obtained  from  the 
following: 

min       S; 


subject  to  P(A)  >  Sj)  <  b2i 


where: 


Xt  =  discrete  random  variable  representing  attrition 

demand  for  item  i  during  a  procurement  leadtime 

b2,  =  parameter  chosen  by  the  decision  maker 
for  item  i 

Using  the  reorder  point  R„  place  an  attrition  buy  when  the  total  number  of  unattrited 
units  (working  and  failed,  installed  and  uninstalled)  reaches  or  drops  below  Rr 
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c.     Step  3 

Use  the  Wilson-Harris  EOQ  model  to  choose  procurement  quantity: 


*  ■  Vt 


where: 

A:  =  administration  and  ordering  cost  for  item  i 
D,,  =  attrition  demand  rate  for  item  i 
/.  =  annual  holding  cost  rate  for  item  i 
C,  =  unit  cost  for  item  i 


So  the  initial  buy  quantity  formula  is  R,  +  0'  =  K,  +  S,  +  Q'. 

This  proposed  repairable  item  inventor}'  control  process  works  as  follows. 
To  start  the  system,  purchase  A'  +  S,  +  O'  units.  As  installed  items  fails,  repair  them 
as  soon  as  possible  (i.e.  do  not  use  batch  repair).  When  attritions  in  the  repair  process 
reduce  the  actual  population  size  to  K,  +  S„  place  an  attrition  buy. 

D.     SUMMARY 

In  this  chapter,  we  developed  a  new  inventory  model  based  on  the  M/M/l/K/K 
queueing  theory  and  the  current  repairables  management  system  of  the  Korean  Air 
Force.  In  the  next  chapter,  we  will  compare  the  current  inventory  model  with  the  pro- 
posed model,  in  terms  of  the  numerical  results  of  each,  using  computer  computations. 
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V.     COMPARING  THE  COMPUTATION  RESULTS  OF  THE  MODEL 

A.  INTRODUCTION 

In  this  chapter,  we  will  discuss  a  specific  result  for  the  two  models,  one  is  a  existing 
model  of  the  Korean  Air  Force,  the  other  is  a  proposed  model.  The  results  show  only 
summary  of  the  computer  computation  and  a  comparison  of  the  results.  The  whole 
computation  results  are  shown  in  Appendix  C. 

The  purpose  of  this  chapter  is  to  provide  a  numerical  analysis  between  the  proposed 
model  and  the  current  Korean  Air  Force  model.  The  data  used  in  the  computer  com- 
putation consists  of  a  selection  of  8  line  items  form  the  data  on  50  line  items  which  was 
obtained  from  the  Korean  Air  Force.  Due  to  the  limited  time  available  for  this  research, 
only  8  items  were  selected  for  use  in  the  comparison  of  the  models.  A  selection  of  high 
demand  and  low  demand  items  were  chosen.  The  complete  data  set  is  described  in  Sec- 
tion C  of  this  chapter. 

B.  MEASURE  OF  EFFECTIVENESS  (MOE) 

The  following  MOEs  are  all  reasonable  methodologies  for  evaluating  repairable  item 
inventory  models  for  the  Korean  Air  Force.  However,  we  will  concentrate  on  one  MOE 
for  model  comparison  such  as  operational  availability. 

1.  Fill  Rate 

It  is  computed  by  taking  the  total  number  of  units  demanded  at  a  base  over  a 
fixed  period  time  and  dividing  that  number  into  the  total  number  of  units  issued  at  the 
time  they  were  are  demanded.  Thus,  fill  rate  is  the  percentage  of  demands  that  are  filled 
immediately. 

2.  Operational  Rate 

Operational  rate  is  the  probability  that,  at  any  given  point  in  time,  there  will  be 
no  stockouts  from  base  supply  (backorders).  Operational  rate  is  computed  by  counting 
up  the  length  of  time  (in  days)  during  a  year  that  no  backorders  existed  and  dividing  this 
number  by  365.  This  gives  us  the  percentage  of  time  during  the  year  that  no  backorders 
were  in  existence. 

Operational  rate  has  an  advantage  over  both  the  fill  rate  and  mean  number  of 
backorders  in  that  it  may  be  directly  related  to  the  supply  system's  effect  on  operations. 
However,  it  has  a  disadvantage  over  both  the  fill  rate  and  mean  backorders  in  that  it  has 
a  rather  bothersome  all-or-nothing  character. 
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3.     Mean  Supply  Response  Time  (MSRT) 

MSRT  is  the  mean  time  it  takes  for  the  supply  system  to  respond  to  the  demand 
for  a  replacement  part  or  component.  It  is  obtained  by  taking  each  stockout  that  is  es- 
tablished during  a  period  of  time,  observing  how  many  days  it  takes  to  satisfy  the 
backorders.  adding  up  all  these  numbers  and  dividing  the  sum  by  the  total  number  of 
demands  during  that  period.    It  is  the  average  time  weighted  unit  short  (TWUS). 

MSRT  is  important  by  itself  as  an  indicator  of  success  of  the  supply  system  in 
meeting  response  time  goals.  As  a  measure,  it  considers  both  the  likelihood  of  satisfying 
demands  from  stock  on  hand  and  the  length  of  the  delay  in  satisfying  demands  when  the 
system  runs  out  of  stock. 

■4.     The  Average  Number  of  NORS  Aircraft 

NORS  stands  for  "not  operationally  ready  because  of  supply".  The  average 
number  of  NORS  aircraft  is  computed  by  counting  for  each  aircraft  the  number  of 
NORS  days  during  the  course  of  a  year,  adding  these  numbers  up  for  each  aircraft,  and 
finally  dividing  by  365.  Considering  that  the  purpose  of  a  spare  parts  supply  system  for 
the  Korean  Air  Force  is  to  maintain  the  operational  readiness  of  aircraft,  the  average 
number  of  NORS  aircraft  would  certainly  seem  to  be  a  reasonable  measure  of  effective- 
ness. A  stockage  model  attempting  to  optimize  with  respect  to  a  NORS  measure  would 
require  more  restrictive  assumptions  than  those  optimizing  with  respect  to  fill  rate,  the 
average  backorders,  MSRT.  or  the  operational  rate.  Above  all,  its  lack  of  "separability" 
is  a  cause  of  major  mathematical  problems. 

Something  to  note  is  that  the  operational  sector  of  the  KAF  has  consistently 
suggested  that  the  Korean  Air  Force  supply  system  should  be  managed  in  operational 
terms,  not  by  supply  terms.    That  is,  their  imminent  question  has  been  "How  many 
spares  are  needed  to  provide  an  operational  readiness  of  X  percent  of  the  aircraft?". 
5.     Constraints 

The  budget  available  for  the  investment  in  spares  is  either  a  constraint  in  all  of 
the  models  or  the  objective  is  to  achieve  specified  performance  at  the  minimum  budget. 
The  essence  of  any  inventory  control  problem  is  the  trade-off  between  cost  and  system 
performance.  As  mentioned  in  the  beginning  of  Chapter  II,  the  budget  constraint  is  the 
major  resource  constraint  for  the  Korean  Air  Force. 
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C.     COMPARING  THE  MODEL  RESULTS 
1.     Data  Characteristic 

The  data  used  for  the  model  comparison  are  8  line  items  of  the  50  line  items  in 
the  sample  data.  The  whole  sample  data  is  shown  in  Appendices  A  and  B.  This  data 
consists  of:  a)  quarterly  demand  data  for  19  quarters;  and  b)  quarterly  attrition  data  for 
19  quarters  for  each  item.  The  demand  data  was  analyzed  using  Minitab.  The  mean, 
standard  deviation,  variance  and  squared  coefficient  of  variation  were  calculated  for  each 
item.  Since  data  concerning  the  failure  process  for  each  installed  item  wasn't  available, 
we  estimated  the  mean  item  failure  rate,  a,  by  computing  mean  demand  over  5  years, 
divided  by  number  of  aircraft.  We  estimated  the  annual  item  repair  rate,  n„  by  dividing 
actual  working  hours  per  year  by  the  average  repair  time  in  hours  for  the  item.  The 
current  average  repair  time  for  each  of  the  8  items  in  the  sample  was  obtained  from  the 
Korean  Air  Force  [Ref.  23].  These  values  are  shown  in  the  table  below.  For  the  pur- 
poses of  this  comparison  we  assumed  that  the  random  variables  associated  with  the 
failure  process  and  the  repair  process  had  exponential  distributions. 


Table  4.     MEAN     REPAIR     TIME: 

This  table  shows  the  mean 
active  repair  times  for  the  8 
sample  items 


Item  number 

Mean  Repair 
Time  (hours) 

1 

2.9 

7 

5.3 

o 

3.9 

4 

2.1 

5 

7.5 

6 

6.1 

7 

1.5 

8 

1.2 

The  actual  data  for  an  item  consists  of  19  quarterly  demand  observations  during 
1984-1988.  Items  were  selected  from  three  different  kinds  of  aircraft,  the  F-4D/E, 
F-5A/B  and  F-16A/B.  Appendix  B  shows  the  Minitab  results  for  all  50  items  in  the 
sample  data  set. 
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The  data  used  in  this  research  also  included  quarterly  attrition  quantities  for  19  quarters. 
However,  for  the  model  comparison,  we  selected  8  items."  The  items  which  were  chosen 
had  a  variety  of  annual  demands. 

2.  Comparing  the  Results 

The  complete  computation  results  for  both  models  are  shown  in  Appendix  C. 
A  summary  of  these  results  is  shown  on  the  following  pages.  For  the  Korean  Air  Force 
model,  the  target  operational  availability  is  90  percent.  In  the  Korean  Air  Force,  this 
is  a  fixed  level  of  availability. 

During  the  computational  steps,  the  actual  operational  availability  allows  an 
increase  from  90  percent  to  93  percent,  however  the  proposed  model  didn't  meet  the 
budget  constraint  for  a  94  percent  operational  availability.  Figure  8  shows  the  difference 
of  K  values;  a  comparison  of  the  current  K  values  of  the  Korean  Air  Force  with  min  K 
from  the  proposed  model.  Figure  9  shows  the  difference  in  safety  levels  for  the  proposed 
model  and  the  Korean  Air  Force  model.  In  the  comparison,  the  attrition  protection 
levels  of  0.98  is  fixed.  Appendix  C  shows  the  results  based  on  several  attrition  protection 
levels  between  0.90  and  0.98.  SLQ  is  a  safety  level  for  the  Korean  Air  Force  model,  min 
S  is  the  attrition  safety  level  for  the  proposed  model.  Figure  10  shows  the  difference  in 
procurement  quantities.  For  the  Korean  Air  Force  model,  RO  is  a  procurement  level, 
optimal  Q  is  for  the  proposed  model.  In  the  actual  difference  for  the  two  models,  the 
proposed  model  has  slightly  lower  procurement  quantities  than  the  Korean  Air  Force 
model. 

3.  Description  of  the  Computer  Program 

For  the  Korean  Air  Force  repairable  item  inventory  model,  we  first  computed 
a  demand  forecast  using  a  trend  value  to  guide  our  choice  of  alpha  value  to  be  used  in 
an  exponential  smoothing  forecast  for  the  next  quarter's  demand.  When  the  trend  value 
fell  between  0.9  and  1.1,  we  chose  an  alpha  value  of  0.2,  otherwise  we  used  an  alpha  of 
0.4.  Then  the  demand  forecast  was  used  to  calculate  the  DDR,  DRP,  RCQ,  OSTQ,  SLQ 
and  RO.  The  computational  procedure  for  this  current  Korean  Air  Force  model  is 
briefly  described  in  the  following  paragraph.  The  complete  computational  results  for  the 
Korean  Air  Force  model  are  shown  in  Figures  20  to  27,  and  the  whole  mathematical 


"    Item  numbers  for  the  8  items  chosen  from  the  50  items  are  10,  16,  20,  27,  30,  33,  37  and 
46. 

8   Attrition  protection  level  is  the  probability  that  attrition  will  not  reduce  the  population  of 
item  i  below  K,  when  an  attrition  reorder  is  placed  when  the  population  drops  to  the  K,  +  S,  level. 
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Figure  8.      Difference  of  Total  Stock  Level 
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Figure  9.      Difference  in  Safety  Levels  Between  the  Two  Models 
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Figure   10.      Difference  of  Optimal  Procurement  Quantity 
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process  for  the  Korean  Air  Force  repairable  item  inventory  model  is  described  in  Chap- 
ter II. 

For  the  DDR.  we  computed  the  total  demand  for  each  quarter,  divided  by  365 
days.  The  DRP  equals  the  RTS  divided  by  the  sum  of  RTS.  XRTS  and  COXD.  The 
RCQ  was  computed  by  multiplying  DDR  times  DRP  times  RCT.  This  RCT  refers  to 
the  repair  cycle  time  allowed  for  depot  level  maintenance.  For  these  computations,  we 
used  120  days  for  RCT  because  this  is  the  maximum  allowed  by  the  Korean  Air  Force. 
The  OSTQ  is  computed  by  multiplying  DDR  times  XDRP  times  OST.  For  OST,  we 
used  365  days  which  is  the  maximum  allowed  by  the  Korean  Air  Force.  The  SLQ  is 
computed  by  finding  the  square  root  of  3  times  RCQ  plus  OSTQ.  The  RO  is  computed 
by  adding  SLQ.  RCQ  and  OSTQ. 

For  the  model  comparison,  we  used  two  system  parameters  (SLQ  and  RO)  and 
the  item's  current  K  value.  Using  the  historical  demand  data  obtain  from  the  Korean 
Air  Force,  and  their  current  repairable  item  inventory  as  described  above,  the  requisition 
objective  (RO)  and  safety  level  quantity  (SLQ)  were  computed  for  each  calendar  quarter 
in  19SS  (note  that  the  equivalent  historical  information  was  not  available  from  the 
Korean  Air  Force).  The  current  K  value  in  Figure  8  was  obtained  by  taking  the  sum 
of  RFI,  XRFI  and  the  currently  installed  population,  as  of  Xovember  1988  (previous 
historical  data  was  unavailable).  The  SLQ  for  model  comparison  in  Figure  9  is  the  av- 
erage safety  level  in  1988.   Also,  the  RO  in  Figure  10  is  the  average  RO  level  in  19SS. 

For  the  proposed  model,  we  first  need  to  choose  bu  for  the  first  constraint  set 
for  each  item.  For  bx.,  one  minus  the  current  target  operational  availability  in  the 
Korean  Air  Force  was  used.  To  obtain  the  constraining  total  budget  value  B,  we  com- 
puted the  sum  of  the  current  Koi  value  for  each  item  times  its  unit  cost  plus  the  holding 

a, 
cost  for  each  item,  over  all  8  items.    The  rho  value  ( — )  for  an  item  was  derived  from 

the  average  failure  rate  (a)  of  that  repairable  item,  divided  by  its  average  repair  rate  (n). 

For  the  model  comparison,  a  value  for  the  budget  constraint  had  to  be  obtained. 
We  estimated  a  value  for  this  variable  by  summing  the  following  products  over  all  items: 
Koi  times  unit  cost  plus  annual  holding  cost  for  item  i  (reminder:  Koi  is  the  number  of 
units  of  item  i  currently  owned  by  the  Korean  Air  Force). 

In  step  1,  the  min  K,  values  are  computed  based  on  two  constraints.  First,  we 
calculate  the  probability  that  .V,  is  greater  than  K,  —  D,  for  an  item.  The  starting  value 
used  for  K,  is  D,  (the  size  of  the  currently  installed  population).  For  the  min  K,  in  the 
computer  program,  we  first  compute  P0  of  the  M/M/l/K/K  queueing  system.  And  then, 
we  compute  P{\,  =  n)  using  P0.  In  the  P(X,  =  n)  computation,  n  equals  K,  —  D,.  If  the 
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constraint  on  P(Nt  >  K,  —  D)  isn't  satisfied  by  choosing  A'  =  D,  then  we  increased  A'  by 
1.  These  steps  are  continued,  until  P(X,  >  A'  -  D)  is  less  than  0.10.  At  this  time,  we 
took  the  current  n  value,  and  then  added  Z),  to  get  the  min  A,  value.  After  obtaining  the 
min  A'  values  from  the  first  constraint  set  for  each  of  the  8  items,  the  budget  constraint 
is  checked.  If  these  K,  values  don't  satisfy  the  budget  constraint,  we  go  back  to  the  first 
constraint  set  and  increase  the  bu  values.  In  this  case,  a  method  for  choosing  an  in- 
creased set  of  bu  's  must  be  chosen.  Two  ways  exist  to  do  this.  The  first  method  is  to 
increase  the  bu  values  simultaneously  for  all  values.  The  other  way  of  solving  this 
problem  is  to  apply  different  operational  availabilities  for  each  item.  In  other  words,  the 
operational  availability  chosen  would  depend  on  the  item  characteristics.  The  first 
method  described  above  (simultaneous  increase  of  every  bu)  was  used  in  this  work.  If 
the  min  A*  satisfies  the  budget  constraint,  this  is  min  A",  used  the  model  comparison,  in 
Figure  S.  As  it  turns  out,  the  first  set  of  min  A,  values  met  the  budget  constraint.  So, 
the  min  A'  values  from  the  first  constraint  set  gave  a  feasible  solution  in  the  budget 
constraint.  Actually,  the  budget  constraint  allowed  up  to  93  percent  operational  avail- 
ability for  all  8  items. 

For  step  2,  we  use  the  Poisson  probability  distribution  to  get  minimum  attrition 
safety  levels.  Suppose  that  risk  b2l  =  0.1.  This  would  mean  that  the  attrition  protection 
level  is  0.90.  Therefore,  1  —  b2,  equals  the  protection  level.  For  the  min  S„  we  first  cal- 
culate probability  of  X,  =  0,  and  next  calculate  probability  X,  =  0  and  so  on.  We  sum 
these  individual  probabilities  up  until  the  cumulative  value  meets  or  exceeds  the  pro- 
tection level  (for  the  model  comparison,  we  use  0.98).  This  is  the  min  S,  value  for  the 
model  comparison  in  Figure  9. 

For  step  3,  the  Wilson-Harris  EOQ  formula  is  used.  For  this  formula,  a  per 
order  ordering  cost  of  10  percent  of  the  unit  cost  was  used.  This  is  the  value  which  the 
Korean  Air  Force  currently  believes  to  be  in  effect  [Ref.  1:  p.  65].  Similarly,  a  holding 
cost  rate  of  10  percent  per  year  of  unit  cost  was  used.  The  resulting  optimal  procure- 
ment quantities  for  the  model  comparison  are  shown  Figure  10.  The  complete  compu- 
tational results  for  the  proposed  model  are  shown  in  Figures  28  to  35. 
4.     Summary  of  the  Results 

The  proposed  model  allows  for  a  93  percent  operational  availability  within  the 
budget  constraint.  This  is  an  improvement  of  3  percent  in  operational  availability  with 
the  same  budget  that  was  used  with  the  current  Korean  Air  Force  model  to  obtain  90 
percent  operational  availability. 
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D.     SUMMARY 

In  this  chapter,  we  compared  two  models  based  on  computer  computations  using  a 
sample  data  set.  The  proposed  model  has  shown  an  improvement  over  the  current 
model.  However,  it  doesn't  mean  the  proposed  model  is  better  for  the  Korean  Air  Force. 
A  tremendous  effort  would  be  required  to  actually  implement  this  model  in  the  Korean 
Air  Force.  Furthermore,  additional  research  should  be  done  to  explore  the  performance 
of  the  proposed  model. 


71 


VI.     SUMMARY  AND  CONCLUSION 

For  this  research,  there  were  three  major  objectives.  The  first  was  to  review  the  lit- 
erature to  understand  the  current  Korean  Air  Force  model  and  other  models,  such  as 
METRIC,  Mod-METRIC,  for  controlling  the  stockage  decisions  in  multi-item,  multi- 
echelon  inventory  systems  involving  repairable  items.  The  second  objective  was  to  de- 
velop a  new  model  for  the  Korean  Air  Force  inventory  management  system.  The  last 
objective  was  to  evaluate  this  new  model  by  comparing  the  current  Korean  Air  Force 
model  with  this  proposed  model. 

For  the  first  objective,  we  studied  the  METRIC  family  model,  ACIM  model  and 
U.S  Navy  model.  We  then  reviewed  their  basic  assumptions  and  mathematical  ap- 
proach. For  the  second  objective,  we  reviewed  the  Korean  Air  Force  inventory  man- 
agement system  concept  and  conditions,  and  developed  a  proposed  model.  In  order  to 
accomplish  the  last  objective,  sample  data  from  the  KAF  was  used  in  comparing  the 
Korean  Air  Force  model  and  proposed  model. 

Analysis  of  the  sample  data  revealed  that  the  proposed  model  yielded  some  im- 
provement in  operational  availability  and  safety  level  with  the  same  budget.  As  men- 
tioned in  Chapter  V,  the  proposed  model  requires  further  research.  Currently  the 
research  is  hard  to  do  because  of  time  constraints  and  data  available  from  the  Korean 
Air  Force.  So,  additional  research  should  be  done  to  eliminate  these  shortcomings  in  the 
proposed  model.  Also,  for  the  implementation,  we  need  more  data  points  to  obtain 
better  distributions  for  the  demand  process  and  repair  time  process. 
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APPENDIX  A.     SAMPLE  DATA  LIST  1 
A.     SAMPLE  DATA  DEMAND  DISTRIBUTION  LIST 

Table  5.     DEMAND  DISTRIBUTION  LIST 
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0 

0 

2 

19S5 

Total 

0 

0 

0 

0 

0 

0 

0 

0 

0 

5 

S6-1 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

86-2 

0 

o 

0 

0 

0 

0 

0 

0 

0 

2 

86-3 

0 

0 

0 

1 

0 

1 

0 

0 

1 

0 

86-4 

0 

0 

0 

0 

0 

1 

0 

0 

2 

1 

3 

1986 
Total 

0 

0 

0 

3 

0 

2 

0 

0 

3 

5 

S7-1 

0 

2 

0 

2 

1 

0 

0 

0 

5 

2 

S7-2 

2 

1 

0 

4 

0 

1 

1 

1 

7 

0 

87-3 

0 

0 

0 

0 

0 

0 

0 

1 

2 

0 

87-4 

1 

0 

1 

5 

0 

2 

0 

1 

2 

4 

1987 
Total 

3 

3 

1 

9 

1 

3 

1 

3 

16 

6 

88-1 

o 

2 

1 

1 

0 

0 

7 

0 

4 

3 

SS-2 

0 

1 

0 

3 

0 

3 

0 

1 

6 

0 

88-3 

1 

0 

0 

5 

1 

0 

0 

1 

4 

2 

198S 
Total 

1 

3 

1 

9 

1 

-> 

j 

2 

2 

14 

5 

73 


Table  6.     DEMAND  DISTRIBUTION  LIST 


QTR 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

84-1 

0 

0 

0 

1 

0 

6 

1 

2 

1 

0 

84-2 

0 

0 

0 

2 

1 

1 

0 

4 

7 

3 

84-3 

0 

0 

0 

2 

0 

0 

0 

0 

2 

0 

84-4 

0 

0 

0 

2 

0 

2 

0 

2 

5 

2 

19S4 
Total 

0 

0 

0 

7 

1 

9 

1 

8 

15 

5 

85-1 

1 

0 

4 

0 

1 

1 

0 

8 

4 

2 

85-2 

0 

0 

0 

5 

0 

1 

0 

0 

3 

0 

85-3 

0 

0 

0 

2 

0 

3 

0 

0 

1 

2 

85-4 

0 

1 

0 

3 

1 

0 

0 

0 

1 

3 

1985 
Total 

1 

1 

4 

10 

2 

5 

0 

8 

9 

7 

86-1 

0 

1 

0 

4 

0 

0 

2 

2 

3 

4 

86-2 

0 

0 

0 

1 

1 

2 

0 

5 

^ 

j 

0 

86-3 

0 

0 

0 

3 

3 

0 

0 

1 

5 

3 

86-4 

0 

0 

0 

2 

1 

1 

0 

1 

1 

2 

1986 
Total 

0 

1 

0 

10 

5 

3 

2 

9 

12 

9 

87-1 

4 

0 

0 

5 

2 

1 

0 

1 

0 

0 

87-2 

1 

1 

0 

4 

5 

0 

0 

2 

4 

0 

87-3 

5 

0 

1 

1 

j 

1 

0 

0 

1 

1 

87-4 

2 

0 

0 

3 

2 

1 

0 

0 

2 

2 

1987 
Total 

12 

1 

1 

13 

12 

3 

0 

3 

7 

3 

8S-1 

1 

0 

1 

4 

3 

0 

1 

0 

2 

0 

88-2 

2 

1 

0 

2 

2 

1 

0 

1 

3 

0 

88-3 

1 

0 

0 

4 

0 

0 

0 

0 

1 

2 

1988 
Total 

4 

1 

1 

10 

5 

1 

1 

1 

6 

2 

74 


Table  7.     DEMAND  DISTRIBUTION  LIST 


QTR 

21 

-)-> 

23 

24 

25 

26 

27 

28 

29 

30 

S4-1 

0 

2 

0 

1 

9 

4 

36 

S4 

4 

4 

S4-2 

1 

0 

0 

1 

5 

5 

44 

66 

7 

0 

84-  3 

1 

0 

0 

0 

11 

4 

52 

43 

1 

0 

84-4 

3 

2 

4 

3 

1 

0 

55 

32 

7 

3 

19S4 
Total 

5 

4 

4 

5 

26 

13 

1S7 

225 

19 

7 

S5-1 

1 

4 

1 

10 

0 

2 

101 

71 

2 

3 

85-2 

1 

1 

0 

0 

10 

7 

67 

59 

0 

0 

85-3 

1 

1 

0 

5 

5 

5 

2S 

50 

5 

0 

85-4 

0 

1 

1 

4 

4 

1 

42 

25 

2 

2 

1985 

Total 

3 

7 

2 

19 

19 

10 

238 

205 

9 

5 

86-1 

2 

5 

4 

5 

7 

6 

61 

40 

0 

3 

86-2 

4 

3 

-^ 

j 

1 

8 

1 

17 

39 

4 

1 

S6-3 

0 

2 

1 

1 

S 

1 

95 

12 

5 

0 

86-4 

1 

4 

0 

4 

4 

3 

47 

73 

0 

5 

19S6 
Total 

7 

14 

8 

11 

27 

11 

220 

164 

9 

9 

87-1 

5 

6 

2 

2 

10 

2 

19 

50 

11 

'  3 

87-2 

4 

1 

0 

2 

7 

7 

24 

37 

7 

0 

87-3 

1 

0 

0 

0 

25 

21 

4 

0 

87-4 

0 

1 

7 

0 

0 

4 

70 

24 

13 

4 

1987 
Total 

10 

11 

9 

4 

20 

13 

138 

132 

35 

7 

8S-1 

0 

5 

1 

5 

9 

10 

42 

31 

5 

2 

S8-2 

6 

2 

0 

6 

9 

1 

21 

18 

9 

0 

88-3 

0 

3 

2 

2 

5 

1 

^ 

18 

23 

3 

4 

1988 
Total 

6 

10 

3 

13 

23 

14 

81 

72 

17 

6 

75 


Table  8.     DEMAND  DISTRIBUTION  LIST 


QTR 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

84-1 

0 

0 

0 

9 

1 

2 

2 

0 

12 

1 

84-2 

5 

0 

1 

7 

7 

1 

3 

0 

21 

7 

84-3 

4 

0 

2 

9 

0 

1 

0 

0 

6 

3 

84-4 

0 

0 

1 

18 

2 

9 

3 

5 

2 

1 

19S4 
Total 

9 

0 

4 

43 

10 

13 

8 

5 

41 

12 

85-1 

12 

0 

1 

17 

10 

4 

2 

1 

17 

9 

85-2 

20 

3 

1 

2 

26 

3 

0 

1 

9 

0 

85-3 

1 

0 

3 

2 

32 

1 

1 

0 

13 

5 

85-4 

14 

0 

0 

7 

25 

3 

2 

2 

5 

1 

1985 
Total 

47 

3 

5 

28 

93 

11 

5 

4 

44 

8 

86-1 

17 

1 

0 

4 

30 

5 

2 

4 

10 

0 

86-2 

25 

0 

1 

2 

11 

0 

5 

0 

0 

13 

S6-3 

0 

1 

0 

9 

17 

0 

5 

0 

5 

7 

86-4 

21 

2 

1 

0 

4 

6 

0 

3 

6 

6 

1986 
Total 

63 

4 

7 

4m 

15 

62 

11 

12 

7 

21 

26 

87-1 

5 

0 

0 

4 

9 

3 

1 

1 

j 

5 

87-2 

5 

1 

2 

5 

5 

3 

2 

1 

2 

5 

S7-3 

0 

0 

2 

5 

2 

1 

0 

0 

7 

5 

87-4 

4 

0 

0 

0 

17 

7 

3 

0 

4 

3 

1987 
Total 

14 

1 

4 

14 

33 

14 

6 

2 

14 

IS 

88-1 

10 

0 

0 

0 

0 

1 

0 

0 

0 

0 

88-2 

3 

0 

1 

0 

10 

1 

0 

0 

4 

3 

88-3 

7 

0 

1 

5 

4 

2 

4 

0 

5 

0 

1988 
Total 

20 

0 

2 

5 

14 

4 

4 

0 

9 

3 

76 


Table  9.     DEMAND  DISTRIBUTION  LIST 


QTR 

41 

42 

43 

44 

45 

46 

47 

4S 

49 

50 

84-1 

0 

0 

1 

4 

-: 

14 

28 

5 

3 

13 

84-2 

0 

0 

0 

9 

1 

0 

30 

6 

3 

9 

84-3 

0 

0 

0 

9 

9 

0 

17 

11 

9 

10 

S4-4 

0 

0 

0 

10 

2 

11 

29 

12 

3 

17 

19S4 
Total 

0 

0 

1 

32 

25 

25 

104 

34 

18 

49 

S5-1 

s 

1 

L. 

5 

1 

5 

23 

8 

5 

15 

85-2 

s 

1 

2 

5 

1 

10 

23 

8 

5 

15 

85-3 

1 

- 

0 

8 

4 

0 

13 

13 

0 

S 

S5-4 

3 

■% 

J 

4 

5 

2 

6 

11 

7 

4 

8 

19S5 
Total 

20 

12 

S 

23 

8 

21 

70 

36 

14 

46 

86-1 

4 

0 

0 

10 

0 

9 

7 

12 

2 

20 

86-2 

7 

4 

1 

13 

11 

0 

19 

10 

3 

0 

S6-3 

1 

1 

1 

7 

9 

0 

18 

5 

7 

6 

86-4 

4 

2 

3 

5 

5 

11 

6 

16 

0 

16 

19S6 
Total 

16 

7 

5 

35 

25 

20 

50 

43 

12 

42 

87-1 

0 

5 

0 

1 

6 

9 

5 

9 

0 

8 

87-2 

0 

1 

1 

1 

6 

7 

9 

10 

6 

1 

87-3 

0 

1 

0 

9 

2 

0 

8 

3 

3 

9 

S7-4 

2 

0 

0 

6 

4 

13 

8 

11 

0 

15 

198" 
Total 

2 

7 

1 

17 

18 

29 

30 

33 

9 

55 

88-1 

0 

0 

0 

0 

5 

0 

7 

11 

0 

4 

88-2 

0 

4 

1 

3 

7 

0 

5 

4 

5 

7 

8S-3 

1 

0 

0 

3 

1 

12 

14 

7 

1 

6 

19S6 
Total 

1 

4 

1 

6 

8 

12 

26 

22 

6 

17 

77 


APPENDIX  B.     SAMPLE  DATA  LIST  2 
A.     SAMPLE  DATA  DISPOSAL  RATE  DISTRIBUTION  LIST 

Table  10.     DISPOSAL  RATE  DISTRIBUTION  LIST 


QTR 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

S4-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

84-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

S4-3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

84-4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

85-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

S5-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

85-3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

85-4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

86-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

86-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

86-3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

86-4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

S7-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

87-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

87-3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

87-4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

88-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

S8-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

88-3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

P(  successful 

repair) 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.92 

78 


Table  11.     DISPOSAL  RATE  DISTRIBUTION  LIST 


QTR 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

84-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

S4-2 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

S4-3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

S4-4 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

85-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

85-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

S5-3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

85-4 

o 

0 

0 

1 

0 

0 

0 

0 

0 

0 

86-1 

o 

0 

0 

1 

0 

0 

0 

0 

0 

0 

86-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

86-3 

0 

o 

0 

0 

0 

0 

0 

0 

1 

0 

86-4 

o 

0 

o 

0 

0 

0 

0 

0 

0 

0 

87-1 

o 

0 

0 

0 

0 

0 

0 

0 

0 

0 

87-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

87-3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

87-4 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

88-1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

8S-2 

0 

o 

0 

0 

0 

0 

0 

0 

1 

0 

S8-3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

P(  successful 
repair) 

1.00 

1.00 

1.00 

0.92 

1.00 

0.94 

1.00 

1.00 

0.94 

1.00 

79 


Table  12.     DISPOSAL  RATE  DISTRIBUTION  LIST 


QTR 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

84-1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

S4-2 

0 

0 

0 

0 

1 

0 

0 

2 

0 

0 

84-3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

84-4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

85-1 

0 

0 

0 

0 

0 

0 

2 

0 

1 

0 

85-2 

0 

0 

0 

3 

0 

0 

0 

1 

0 

0 

85-3 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

85-4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

86-1 

0 

0 

0 

4 

0 

0 

0 

0 

0 

0 

86-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

S6-3 

0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

86-4 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

87-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

87-2 

2 

0 

0 

0 

1 

0 

0 

1 

0 

0 

87-3 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

87-4 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

88-1 

0 

0 

0 

1 

0 

1 

1 

1 

0 

0 

88-2 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

SS-3 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

P(successfu 
repair) 

0.94 

1.00 

0.96 

0.97 

0.97 

0.97 

0.99 

0.99 

0.96 

1.00 

80 


Table   13.     DISPOSAL  RATE  DISTRIBUTION  LIST 


QTR 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

S4-1 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

84-2 

o 

0 

0 

0 

0 

0 

0 

0 

2 

0 

S4-3 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

S4-4 

0 

o 

0 

0 

0 

0 

0 

0 

0 

0 

S5-1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

85-2 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

S5-3 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

85-4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

S6-1 

0 

0 

0 

0 

1 

0 

1 

0 

0 

0 

86-2 

1 

o 

0 

0 

0 

0 

1 

0 

0 

") 

86-3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

S6-4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

S7-1 

0 

o 

0 

0 

0 

0 

0 

0 

0 

0 

8  "-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

8  7-3 

0 

0 

0 

o 

0 

0 

0 

0 

1 

1 

87-4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

88-1 

1 

0 

0 

0 

0 

0 

0 

0 

o 

0 

8S-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

SS-3 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

P(successfu 

repair) 

0.98 

1.00 

0.99 

1.00 

0.99 

1.00 

0.99 

1.00 

0.96 

0.94 

Table  14.     DISPOSAL  RATE  DISTRIBUTION  LIST 


QTR 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

84-1 

0 

0 

0 

0 

0 

2 

0 

0 

0 

1 

84-2 

0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

84-3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

84-4 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

85-1 

0 

0 

0 

0 

0 

1 

0 

0 

1 

0 

85-2 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

85-3 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

85-4 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

86-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

86-2 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

86-3 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

86-4 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

87-1 

0 

0 

0 

0 

0 

1 

0 

2 

0 

0 

87-2 

0 

0 

0 

0 

0 

1 

0 

0 

2 

0 

87-3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

87-4 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

88-1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

88-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

8S-3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

P(successfu 

repair) 

1.00 

0.97 

1.00 

0.97 

1.00 

0.99 

0.99 

0.96 

0.92 

0.98 
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B.     SAMPLE  DATA  ANALYSIS 


Table   15.     DEMAND  DISTRIBUTION  ANALYSIS  LIST 


No 

Mean 

Mean: 

S.D 

Var 

Q 

1 

0.364 

0.133 

0.674 

0.454 

3.43 

2 

0.545 

0.297 

0.820 

0.672 

2.26 

3 

0.182 

0.033 

0.405 

0.164 

4.97 

4 

1.909 

3.644 

1.700 

2.S90 

0.79 

5 

1.909 

3.644 

1.700 

2.890 

0.79 

6 

0.727 

0.529 

1.009 

1.018 

1.92 

- 

0.2  "3 

0.075 

0.647 

0.419 

5.59 

8 

0.455 

0.207 

0.522 

0.275 

1.32 

9 

3.091 

9.554 

2.343 

5.4S9 

0.57 

10 

1.316 

1732 

1.293 

2.940 

0.98 

11 

0.S95 

0.801 

1.449 

2.099 

2.62 

12 

0.211 

0.044 

0.419 

0.176 

3.99 

13 

0.316 

0.099 

0.946 

0.895 

8.97 

14 

2.632 

6.927 

1.422 

2.022 

0.29 

15 

1.316 

1732 

1.176 

2.055 

1.19 

16 

0.S42 

0.709 

0.834 

0.696 

0.98 

17 

0.211 

0.045 

0.535 

0.286 

6.36 

IS                  1.526 

2.329 

2.118 

4.486 

1.93 

19                 2.579 

6.651 

1.805 

3.258 

0.49 

2"                 1.367 

1.S69 

1.342 

1.837 

0.9S 

21                  1.632 

2.663 

1.S62 

3.467 

1.30 

22 

2.421 

5.S61 

1.742 

3.035 

0.52 

j-i 

1.444 

2.0S5 

1.947 

3.791 

1.82 

24 

2.737 

7.491 

2.642 

6.967 

0.93 

25 

5.362 

31.719 

3.609 

13.025 

0.41 
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Table  16.     DEMAND  DISTRIBUTION  ANALYSIS  LIST 


No 

.Mean 

Mean2 

S.D 

Var 

*-* 

26 

3.211 

10.311 

2.594 

6.279 

0.65 

27 

45.470 

2067.521 

24.880 

619.014 

0.30 

28 

42.000 

1764.000 

20.700 

428.490 

0.24 

29 

4.684 

21.939 

3.667 

13.447 

0.61 

30 

1.789 

3.201 

1.782 

3.168 

0.99 

31 

8.050 

64.810 

7.920 

62.726 

0.97 

32 

0.421 

0.177 

0.838 

0.702 

5.65 

0.895 

0.801 

0.875 

0.766 

0.96 

34 

5.530 

30.581 

5.230 

27.353 

0.S9 

35 

11.160 

124.550 

10.420 

108.576 

0.87 

36 

2.789 

7.779 

2.463 

6.066 

0.78 

37 

1.526 

2.329 

1.467 

2.236 

0.96 

38 

0.947 

0.S97 

1.508 

2.274 

2.54 

39 

6.790 

46.104 

5.700 

32.490 

0.70 

40 

3.6S4 

13.572 

3.284 

10.785 

0.79 

41 

2.053 

4.215 

2.838 

8.054 

1.91 

42 

1.579 

2.493 

2.063 

4.256 

1.71 

43 

0.842 

0.709 

1.167 

1.362 

1.92 

44 

6.053 

36.639 

3.353 

12.496 

0.34 

45 

4.632 

21.455 

3.077 

9.468 

0.44 

46 

5.630 

31.697 

5.380 

30.112 

0.95 

47 

14.890 

221.712 

8.720 

76.038 

0.34 

48 

8.842 

7S.181 

3.404 

11.587 

0.15 

49 

3.105 

9.6414 

2.622 

6.8759 

0.71 

50 

9.420 

88.736 

5.370 

28.837 

0.35 
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APPENDIX  C.     COMPUTER  PROGRAMS  AND  RESULTS 


program    ROKAF    (Input. output): 

typ-s 

tnblO    =    nrrny    (1..19.1..0)    of    real; 

tob20    =    nrroy    [I..Q]    of    ro»l; 

tob3Q    =    army    l  1  .  .  0  ,  1  ..  19  ,  1  .  .  0  ]    of    real; 

vo  r 

c lik    !    boo  Iron; 

l.J.k     :     integer; 

in_lnb     :     tnblO; 

Dnil_tnb     :     toblO; 

DUP_tnl>    :     tub  10; 

nCQ.t.ob    :     tnblO; 

0STO_tnb     :     tnblO; 

SI.O_t«b     :     toblO; 

HO_tnb       :     tabid; 

Al.l'_tab     :     tnb?0; 

mil!)    tnb 

tabZfl; 

FCST_tnb 

tnb-.MI; 

EXP_Leb 

tnb 10: 

DIS_tnb 

tnblO; 

TEIIP_tnb 

tnb3l); 

PtlOCF.DUnE      titledi:  integer); 

l)C(J  In 

enso    k    oT 

1          ii  r  1  to  1  ii    ( 

****      TRENfJ      Alt  AI.YSI  S         **  +  ** 

) 

2     :     ur  i  to  In    ( 

»***       EXPONENTIAL         SHOO  Til  I  MG 

****♦•   )  ; 

3         w  r  1  i.o  1  ii    ( 

»rn                         |)DI1                      ***** 

) 

'1         ii r  1  Lc  In    < 

****                         D  Ht*                      ***** 

> 

5    :    Mr  1 teln    ( 

mt                         HCO                      ***** 

) 

6     !     "rite  In    ( 

tiit                         OSTQ                    ***** 

) 

7        ur i  te in    ( 

****                         Sl.Q                      ***** 

) 

0    :    write  In   <"    ****                       HO                      ***** 

) 

end    (enso    It    ) 

and  ; 

PROCEDURE    rec    ( l.j ,k: integer ;var   TEIIP_tob    :  tnb30    ); 

ber,  i" 

c  n  s  o    k    o  1" 

1    :    TEllP_tnb[k.  i  .  j  ] 

=    IH_tnb    [i.J]; 

2     :    TENP_tab[k.  1  ,  i  ] 

=    EXP_tab    [l.j] 

3    :    TEIIP_tal>lk.  i,  j] 

=    UUH    tnb    [i  ,  j) 

4    :    Ti:iir_tob[lt,  l,.i  1 

=    DRP_tnb    [i.J] 

5     :    TI!ir_tol)[k,  i,  j  ) 

=    liai_tab    [  i.j  ] 

6     :    TENP_tah[k,  l.j] 

=   OST«_tnb[i, j] 

7    :    TEIIP_tab[k,  I  ..i] 

=    Sl.q_tab    ti.  j  J 

8    :    TEIIP_t«b[k,i.j] 

=    UO_tnb    [i.jj: 

end     (case    k) 

end  ; 

'igure   1 1.      Computer  programs  for  the  KAF  Model 
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(**•«**      MAIM         IMtOGKAII      ****) 

begin 

(*>>>>       DATA    lliniT    ROUTINE         <<<<*) 

ur  itcln    (  -           DF.HAIII)      DATA    !  !  !       '  )  ; 

for    1     :=    1    to    19    do 

b  c  p,  i  n                                                                                     :! 

ur  i  Lrln    (  '     FOR       D'  .  i  )  ; 

for    .1     :=    1     to    0    r|o 

rend    (  in_tab(  i , .)  ) ) ; 

end  ; 

writcln    ('       DISPOSAL      RATE       !!!     '); 

for    i     : =    1    to    19    do 

bcf»  in 

wri  teln    (  "     FOU       UK"  ,  i); 

for    j     : -     1    to    0    do 

rend    (DlS_tnb    [     i  .  j  1  )  ; 

ond  ; 

(*>>>>      TKKHD         ANALYSIS    <<<<*) 

borj  i  ii 

lor    i     :=    1    to    ()    do 

bef.  i  n 

TntlD_tnb    [il     :=    (  (  (Jll.l.nbf  J().  i]  Hll_lab|  19.  i  J  )»2)/(  Ul_tnbf  16, 

i  ]  i  lll_tnb|  17, 

U* 

lH_tnb[ 10, i 1 

♦    III.. tub [  10.il)): 

if    (THUD..  Lnh    |i|     >=    0.9)       mid    (THUD    tnli[ij     <=    1.1)    then 

ALP_tab    [j|     :=    U.Z 

cine 

ALP_tnli[  i  1     :=    O.'l  : 

end  ; 

oml ; 

(*>>>>         b'XrOIIKUTl  Al.      SMOOTH  1.11(3            <<<<*) 

for     i     :=    1    to    I J    do 

KXP_tnb    (1  .  i  )    :=    HL.tnb    f  1 .  i  J  ; 

for    i     :=    2    to    19   do 

beg  in 

for    ,i     :=     1     to    0    do 

EXP_tnb    [i,.i]     :=    (ALP    tnb    [.j]    *    lll..tab    [  i  ,  j  ]  )    +    ((     1    -    ALP 

_tnb    [j]) 

*    L'XI'_.tnb    (i-l.jj); 

ond  ; 

Tor    i    :=    1    to   0   do 

FCST_tob    [i]     :=    (ALP. tnb    [ij    *    IH^tnb      [19. i])    +    (d    -    ALP 

_tnb    [I]) 

1    F.XFJ.nb    llS.i])'; 

(*>>>>         DAILY      DF.IIAIID      RATI-            <<<<*) 

for    1    :=    1       to    19      do 

begin 

forj     :=    1    to    0    do 

DDR_tab    [    i.JJ    :=    IH_tab    [    i,.i]    /    91; 

end  ; 
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(«>>>>  DEPOT    REPAIR    PERCENT  <<<<*) 

for     I      :=     1     to     19    do 
be*;  in 

for    .i      :=     1     to    8    do 

iT    OIS_tab    [i.j]     =    U.O  tl.cn    DItP_tab    [i.j]     :=     1.0 

elr» 
be?  i  M 

iT     ll.'_tnb    [i.j]     =    U     tlirn    DRP_tab    [i.j]     :=    0.0 
else 

URP_lnb    (i.j    ]     :=    l>IS_tab    [i.j]    /    lH_tob    [i.j]; 

end  ; 
ir    (IHS_tab    [i.jj    =    1)    and    (i::_lab    [i.j]    =    D     tlien    DRF.tnb    [i.j]     :=    0.0 
i ; 

(♦>>>>  KCQ    A    (ISTH    SLU    l<0  c<<<*) 

f t     i      :  =     I     to     10    do 
l)T  i  ii 

for    j : =     1     to    0    do 

beg  in 

RCQ_tab    [i.j]     :=    DOR.tafa    [i.j]    »    URF_tnb    [i.j]    *     120; 

USTC|_tab    [i.ij     :=    M)H_tnb    [i.j]    *    (     1    -    DRP_tab    [i.j])    *    305; 

SI.Q_l«b    I  i.jj     :=    r'iri.    (     3    *     (RCCJ„tab    [i.j]     +    OST0_tnb     [i.j])); 

R0_tab    [i.jl     :=    SHl_tnb[i,j]     ►    RCQ_t ab [ i . j ]    »    USTQ_tab    [i.jj: 

end  ; 
end  ; 

(*>>>>  Pn  HIT    MJWT1RE  <<<<►) 

clik     : =     Tn '      - 

Tor    k    :=    1    lo   n  ,jo 
b>- 

ur  i  I  elii; 

v  r  i  tr  I  ii  ; 

write  In; 

tilln    (k): 

write lii    (  •    UTIl  I  2  :) 

write  In    (' 

'  ); 

Tor     i     :=     1     to     10    tta 
bcfl  in 

write   (  "      II". i); 

f nr     i      : =     1     to    4    da 

bi'i!  iii 

---     (i.j  .k.TKIIP_lr 
write    (  "     "  .TF.Iir_tnb[k,  i.j]); 
end ;  i 

write  In; 
pud  ; 
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if     li    -     I       then 

be;*  in 

wr  i  to  In ; 

wr  i  teln; 

write    ( '       ALP" ); 

Tor    j     :=     1    to    '1    do 

wri te    ( '     ' ,ALP_tab    [j] ); 
wr  i  teln  ; 
end  ; 

if    k    =    2       Llion 

bert  in 

wr  i te  In  ; 

wr  i  te  In  ; 

write    (  '       FCT'  ); 

for    j     :=     1    to    4    do 

write    ('     ' ,FCST_tab[j ] ); 

wr  i  te ) n ; 
end  ; 

wr  i  tn  In  ; 

w  r  i  t  e  1 1 1 ; 

writoln    (  '    QTR 

5                                               6 

7 

0); 

lor    i    :=    1    to    19   do 

bcrtin 

wr  i  to    (  '    D'  .  i ) ; 

Tor    .i     :=    T«    to    0    do 

hei!  j  ii 

roc    (  i  ..j  ,l(  .TKMP  .tnb) 

wiUo(  -     •  ,  TKMP..  tnb  1  It 

i  .  .i  1  )  : 

cud  ; 

wr  i  te  In  : 

pud  ; 

ir       It    =     1          Uion 

bor!  i  n 

wr  i  to  1  ii  ; 

wr  i  to  In ; 

wr  tl.e(  '    At.P'  ) : 

f n r    .i     :  =    S    to    0    do 

write(  '     '  ,AI.P_l.nb[.i  J  )  ; 

writoln; 
end  ; 

iT    It    =    2    thou 

bog  in 

wr  i  te  In  ; 
i            wr  i  te  J  ii  ; 

write(  '    FCT'  ); 
for   j    :=   5   to      0   do 
write('     '  ,FCr>T_tobU]); 
writoln; 
end  ; 

end  ; 

end  . 

Figure  14.      Computer  programs  for  the  KAF  Model 
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program   test; 
typi 

nrray_a  =  array   [J.. 8]   of  real; 

(»rray_K   =  array   T1..0]   of    intorjer; 

orray_x  =   array    [1..9J   oT    reel; 
var 

nrr.n    :    arrny_n; 

arr_.K    :    array.K; 

nrr_Kc:    array_K; 

arr.D    :    array_K; 

nrr.rho    :    array.n; 

nrr_DL:    array_a; 

nrrj<    :    array.. x; 

nrr_l)n    :    nrray_a: 

nrr  _Q      :    ariay_n; 

orr_S    :    array. K; 

i.total .0 ,n. itcm_no,x.q_tcmp   :    Integer; 

("act. sum, PU.exo. cord., prot_level    :    real; 

do_a»'.aiti    :    boolean; 

conr;l. 

BiKlC.ct   =    lDO-lOOOU: 
IC  =  0.1; 
0/\  =  0.9; 

(*>»>     STEP  I     <<<<*) 

function  rtet  JU(Kc;  interter; rho: real): real ; 
v»i  r 

k    :    inl.o'.or; 

\ 
bertin 

I  net   :=  1: 
pi  mm    '.-  0; 

lor   l(    :=  tj   to  Kc  iln 
bcrt  in 

for    i:-   Kc  down to  Ko-k»l  do 
fact    :=    C.not  *    j *  rho : 
(   uriteJu('k=    '.I'..'     fact=    '.hid.);) 
rnni:-   Mini    I    (  o<; t. : 
{    write  Jn(  'r.mn   =    '  ,  rum  ) ;  ) 
rani.    :=    J  ; 
niv I ; 
rtctj'0    :=    l/nufli; 
writeln(  TO  =    ' ,1/rajm); 
curl; 
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procedure  initinl_n(vor  nrr_n:nrmy..n) ; 
beft  in 

nrr..n(lj  :=    l]670.03rl.  1; 

nrr.nf?.)  :r  3G70.nO* I.  1; 

nrr_n|  3)  :=   5M0'1 . (JO ^  1 .  1; 

nrr..n|<l]  :=    1  KJ01  .  0?.*  I .  I; 

nrr..n|r>]  :=   32-110.00*1.  1; 

nrr_nlC|  :=  70G7.05U.1; 

orr_r>!7]  :=    110C0.72»1.  1; 

nrr.nlO]  :=   51013. 7H  1 .  I ; 
end; 


procrdnro  dirplny_minK(a 
bnf.  in 


urii.c)n( 
mile  t  ii  C 
wr  i In  ln( 
ur i In |n( 
wr  i  l.cln( 
ur  i  l.c  ln( 
wr  i  LcJn( 
write ln( 
end ; 


miii  K 1  = 

min  Y.2  = 

miii  K3  = 

min  I'.-l  = 

min  K5  = 

min  KH  = 

min  R7  = 

min  KO  = 


r_K    :    arr;.iy_K); 


.nrr.Kri]) 
.arr.KlZl) 

.nrr.K|3)) 
.nrr.KI-ll) 
.»rr_K|5J) 

,nrr_KlGl) 
.nrr_K[7)) 
.nrr_K[0J) 


function   n I.I i  power (L'l..    :    renl'.powrr    :    integer ) :  r.eol ; 
bertin 

iT   pounr   =  M   then 
iith.powci'    :=    1 
olee 

nl.li. power    :  =   l»lj  *   nLli  poucr(DI.,, power- 1  ) ; 
riul; 


I'linctirti   ["{»:( fium    :    inl.c/ier):  inl.cf.er; 
bef*  in 

iT   mini    -  IJ    llii-n 

lac    :=    I 
i*  Ire 

I'lic    :=   Hum    '     I  noOnim-]  ) ; 
end; 


fm  ic  Li  on  jiet.od'l,:  renl  ;proU.Jevel :  re.-il ):  intci'rji* ; 
vnr 

I'.r.um    :    ten  1 ; 
i    :    inler'.er; 

bcr!  in 

t'.r.nm    :=   O.U; 

i    :=  0; 

rercnt  < 

P..cnni    :=  r'.r.nm    ►  nt.li_p<>u?r(DI»,  i  )•  (cr.p(  -0L.))/fac(  i ) ; 
i    :-    i    I     1; 
unl.il  P.sum  >:  pro i.   love] ; 
«ot_S   :=   i-l; 
end; 


Figure   16. 
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procedure   i«ct_UHwor   nir.DL    :    niray_a); 

begin 

orr_m.ll] 

=  0.124: 

orr.Dl.m 

=  0.052; 

nrr..DI.[3J 

=  0.15-1; 

orr_I)l-l'Jl 

=  o.30«l: 

nrr.DI.15] 

=  0.230; 

orrJ)l.l6| 

=  0.0rj6; 

■rrJ)!.l7J 

:=  0.310; 

orr.ULlHJ 

:=  U.010; 

end; 

rrocnduro  dir;plny..S(  nrr_S    :    nrrDy_K); 

bcf.  in 

uri  to( '      ' ,nrr_S| 1] ) ; 

writeC      ',nrr_S[2|>; 

uiiL^C      '.arr_G(3]>; 

wriLcC      '.orr.SMJ); 

uriteC       ",arr_S|5J); 

writri('      ',nrr_S|G|)i 

wril<?('       ',nrr_5|7]  >; 

writelnC      \nrr_S(.01); 

end ; 

procedure  gcl_inR<var  nrr_.Kc.an_U   :   array_K;   var  oir_rho   :   array_n); 

bed  in 

arr   Kc|  11    : 

=  00; 

»rr_1lc_2|    : 

-    110: 

orr_Kc[3] 

-    in.';; 

nrr_Kcf41 

=   103; 

«rr.Kc|.r>1 

-   HO; 

nrr.R-MGI 

=    19'); 

nir_Kcf7| 

:=   191; 

nrr„KcL0J 

=   120; 

nrr.HI  1]    :  =   00    ; 

nil'  nm    :=  W    ; 

orr.nm    :-   175; 

prr  l)|  -1|    :=    175: 

orr.lH5|    :="UU   i 

nrr  1)10  |    :;    170: 

nir_IJf7]    :=   175; 

nrrJUOJ    :=  00   ; 

nrr.rliof  11 

=  O.UOGOf) 

nrr_.rliol2] 

=  U. 01 010 

nrr_rho[3J 

=  0.1 1000 1 

orr_rliof '1 1 

=  u.oonvo 

nrr..rhof  5  | 

-  0.0001  ".'1 

nrr_rliolHJ 

=  o.oooo;: 

orr_rlm[7J 

:=  n. 111529 

nrr_rho|01 

=  0.00071 
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nrr_xl 1 J 

= 

U 

"JU 

arr_x[Zl 

= 

U 

91 

nrr_x[3 ] 

= 

0 

92 

nrr_x[<1  ] 

r 

0 

93 

t»rr_*[5] 

= 

0 

9-1 

nrr_x[G] 

= 

0 

95 

nrr_x(7] 

r 

0 

96 

arr_x[01 

=: 

0 

97 

nrr_x[9J 

= 

0 

90 

end; 


begin 

ini tial_a(nrr_n ) ; 
Eet_.init(nrr_Kc.nrr..l),arr_rlio); 

repeat 

do_nr?nin    :=    ["nice; 
item_np    :=   0; 
cxp_cost    :=  0.0; 

repeat 

item. .no    :=    il.rm.no  +    3; 

FU   :=  get_PU(nrr_Kc[ itcmjioj ,arr_rho[ item.noj ); 

Tact    :=    t; 

sum    :=   CI; 

ti    :=  0; 

repeat 

for    i:=   arr_/>cf  il.em_no]   downto   arr_Kc[ item_no]-n»^  do 
fact    :=    fact  *    i*urr_ilio[  item  .no  ] ; 

[net    :=   Cnct^FO;  tnnccr.sary   to  comment) 

write(  '        ii  =    '  ,n ); 

sum:=  sum  +   fact; 

writeln(  '  sum  =    ' ,nnm); 

fact    :=    1; 

n    :=  ii    i-    I; 
nnl.i  1(1-  ruin)    <-   (]    -  ()A); 
write ln( 'min  K  =    '  ,orr.  I)|. item. no  Jiri-1): 
prr.KI  il.rni.noj    :=  err.OT il-rrii  nn|   +  n   -I; 

cap  ..cost   :=  rjxp_cor;t  +  arr_n|.item_no  | :l(ari\_l)|  itcmjio  |  m-1 ); 
until   item_no  =  0; 
if  exp_cont   <=  LWidfJol'.   then 

di?jplny_minK(3rr_K) 
•    olr;9 

tlo_Pflnin    :=    true;  , 

until  not  do_again; 


{+>>>>        STEP   11        <<<.<!) 

wriheln("   prot_level 

writeln( ' 

for  x    :=   1    to  9  do 
beg  i  n 


SI  S2  S3  SI  Sli  S6  S7  SO'); 
-); 


Figure  18.      Computer  programs  for  the  proposed  Model 
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«rite( '      ' ,arr_x[.xj ); 
item_no    :=    0; 
gct_DL(arr_DL); 
repeat 

itcm.no    :=    ilcm_no   *     1; 

nrr_S[ itcm.no]    :=   f.et_S(nrr  J)Ll  itcni_nr.]  ,  orrj([x]  )  ; 
until    itcm.no   =   0; 
d  lru>loy_S(orr_S) ; 
end; 


(*>>>>       STH?  Ill    <«<*) 


nrr_Dnll]    : 

=  0.2: 

nrr_Ua[2J 

=  0.2; 

nrr_Un(3] 

=  0.2; 

nrrJM'I] 

=  0.4; 

Brr_Dnf5] 

=  0.2; 

nrr.DnlGJ 

=  0.2; 

arr_l)ol71 

=  0.6: 

orr.DnlOJ 

:=  0.6; 

urilcln; 

' 

write  J 1 1 ; 

wr i trln    (  ' 

*>>>>    optimal  tntucuwiiirarr  quantity    <<<<+ 

); 

Tor   x   :s  1 

to  U  do 

begin 

nrr.U    [x 

]    :=  aqrt   ((2*   (0.1  *   (arr_pLx]/l.l))*  air. 

.Da 

wri  to    (  ' 

4'  ,x); 

uri tcln( 

'    =    ',nrr_«lxj); 

m  id ; 

uril.eln; 

wr  i  t  c In ; 

wri  l.e  In    (  ' 

K                 S              «               TOTAL 

'); 

up  Until  (  ' 

—             —           —               

'): 

for    x    :=    1 

to  0  do 

bei!  in 

ir  nrrj 

|xj    <=   2    LIk.1I  0    Ifn.p    :-?: 

if  mr.c 

|x|    <  =    I    llieii  q.trMnp    :  =  l: 

botn]    := 

mi.  K[x|   ♦   arr...'?|xl   +  q_l.cmp; 

wr ile(  ' 

\x); 

ir  x  -  ] 

then  writc( '    ' ); 

wr  il.c(  ' 

'.orrji(xj): 

wri tc(  ' 

'.arr_SlxJ): 

writc(  ' 

'  ,q.  temp) ; 

ir  x  =  i 

tlien  Hril.e("    ' ) ; 

wr ilc(  ' 

'.total); 

wri  he In; 

rnd  ; 

end. 

Figure 


19.      Computer  programs  for  the  proposed  Model 
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Figure  26.      Computational  Results  for  the  KAF  Model 


100 


»  »  »  « 

(ITU 

Dl 

02 

03 

01 

05 

HO 

07 

03 

09 

DIM 

01  l 

DI2 

U)  J 

0  M 

D15 

01(3 

1)17 

DID 

UJ3 


no 


t  t  •  «  * 


1 


3.  3O1372G208E 

o . ooouoouoonn 

0  . 0000U00U0UK 
2  .706447907  IK 
2.706447907IE 

o.oooononunoF. 

U.OOOO'JUUOOOK 

2.06202oui29E 

2.  706-147907  IE 

i .321004  lean 

0.0000000000 
3.864S3S24  Tti 

2.706  447007  1 
U.ODOtlllljDDUH 

U  .0000000000 

J  .  3 
2.0G2ll20(i  IT. 

u.ouuuumiuiji! 
6  .-luoe/ouu/u 


i  00 
»00 
•  HO 

•  oo 

•  (JO 

•  00 

•  uu 

•  00 

K  ♦  UU 
IE  1 00 
E«-00 

E»ou 
r»uo 

K»UU 

E  i-OO 

i  HO 

mo 

E+00 


2  .062023U12 
2.70G-1I7907 

u.uoooooooo 

2  .  7U6-1-1  V'J(J7 

2.0G2020012 

0.000000000 

1. 32  109  I  109 

2. 002020012 

2.0G2020012 

0.00000000 

3.  OR'!  53524 

U. 00000000 

1  .32  11)9-1  IG 
2.06202001 

o .  ooooonot 

2  .  70P44790 
0  000000(10 
2.70G44790 

O.UOOOOOOO 


9E»0U 
1EHJ0 
OEmO 
lEfUO 
9K»00 
OK  +  00 
6E+0U 
9E»00 
9E  +  00 
OOEfOO 
70 E+00 
OOK •00 
OOEfOO 
29t>0U 
OOEtOO 
IE»UU 
00E»O0 
7 LEtOO 
out:  too 


UOOOOUUOOOE 

OOOOOOOOUOE 

0620280129E 

OC2U200129E 

32109-1109GE 

3253395206E 

OUOOOOOOOOE 

OB2020U129E 

OOOOOOOOUOE 

0020200129 

0000000000 

70G-1 -179071 

0620200129 

0000000000 

OdOOOOOOOO 

3013226250 

706-1-17907  1 

ODOOOOOOOO 


2.0620280129 


♦  00 

♦  00 

♦  00 

♦  00 
+  00 

♦  00 
+  00 

♦  00 

♦  UU 
K  +  00 
E  +  00 
E  +  00 

r.+  oo 

ft  00 
E»0U 
KiOO 
Et-00 
Ef-UO 
E  +  00 


O.OUOOOOOOUOEmO 
4  .371920300') li»00 
0  .OOOOOOOOOOE'OO 
2  .O62020OI29E  +  OO 
2.06202GU129EmO 
O.0000000UOOE»OO 
2.O620200I29K+OO 
2.7UG4479071EM)0 
6.400070U97GE+OU 

o.ooooooooour.  »oo 

2.70644790/1E+U0 
2.0620280129E»0U 

o.uuuuuuouuoi;mo 
o.ooooooooum:  mo 
1  .32109-1  I896E»00 
2.0620200129E  I  Oil 
0 .UO0000OO00E»0G 
0.  OOOOOOOOUOE  mo 
2.O62O20O129E»OO 


(itn 


01 

0. OOOOOOOOUOE  KIU 

1 

U2 

1 .321894 IQUCEiOO 

1 

nj 

2.  06202UU  129E»0ll 

0 

04 

1  .  ::;:iP'i  i  IOOOEi  no 

'.' 

Ob 

l  .3Z1U94  loauEmu 

1 

DO 

1  .32109-1  IBOOEtUIJ 

1 

07 

2.700447907  lE'OO 

'1 

UP 

O.OOUOUiJMUOOE-tOU 

0 

1)9 

O.UUOOUUUUOi>rmO 

u 

DIU 

1  .  321094  lOUOEtOd 

01  I 

O.UOOOOOOOOOElOO 

D12 

1 .3210fJ410U»:RHIU 

013 

0.OOOOOUOOUOE«00 

DM 

3.3:::.:i-Tj5206EU'O 

D15 

7  .0020280 1 29E  Mill 

DIG 

U.O0OOOOUU0UE*00 

P  17 

0 .  0OOO"OO"J0OE  •  00 

D10 

1.32189  I  10961  iOO 

019 

1.32  1091  lU'JOEmO 

1 

.32189-U036E 
.321894 IP9GK 
.00O0000000E 
.0020X001  :"JE 
.321804 109HE 
.321804 1896E 
.971920:il»BBK 

.  ooooorjnoooE 

.oouooooooor. 

2.0620200129 

o.uooouooooo 

1  .32109-1  1090 
1.321004  1 090 

0.0000000(100 
1  .321094  I'j'JO 

i  .r.:  109-1 1»90 

O.OOUOUOUUUO 
1  .3210941096 

0.0000000000 


♦  00 

♦  00 
i  00 
100 

♦  00 

♦  00 

♦  00 

♦  00 

♦  00 
EmO 
Ei  0(1 
E+00 
KtOQ 
E+00 
KHIU 
EMJU 
K'OO 
t'OO 
E<OU 


1  .32  1094  1B95t>00 
1  . 321894  18UGK«00 
O.OOOOOOOOOOKHIO 
2.706-1 4  790/ lEi  00 
B.  421)2  147684EMJ0 
O.POOOOUUUOOEMJU 
7.9324  10Z390EM)0 
3.3013226Z50Et0U 
3.Q64535247UE+00 
1  .32109-1  IBOGEmo 
1  .321P94  10UGE»00 
3.301322G258E+00 
3.325:)3952Ut5EiOO 
2. 0620200 120E* 00 
O.OOOOOOOOOOEtOl) 
O.OUOOOUOOOUI->00 
3.0G45357.170K+00 
4  .405257302OEi0U 
2.0G20200129E+00 


8.310G009I97K 

o.ooooooooour. 

OOOOOOOOUOE 
5BU5JU1177E 
B645352470E 
4202147604E 

OOOOOOOOOOE 
4052573020E 
9302853O85E 

OdOdnUUOUO 
UUOOOOOOOO 
5090301177 
332739,1173 
920937  100(1 
OOOOOOIJOOO 
0547431634 
0000000000 
0.0000000000 
7.3054905039 


0 

I , 

3. 
6. 

a 

4 
5. 

0. 
0. 

1 . 
1. 

4  . 
0. 
7  . 
0. 


mo 

•  00 

•  00 
I  01 

•00 

•  00 

mo 

'(JO 

'00 

F.HJO 

fi  01] 

F'Ol 

I; '  ()  I 

E«0(1 

EtOO 

E»00 

Ei  on 

I:  •  00 

EK10 


Figure  27.      Computational  Results  for  the  KAF  Model 
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Figure  28.      Computational  Results  for  the  proposed  Model  (OA=  0.90) 
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Figure  30.      Comj5iitation.il  Results  for  the  proposed  Model  (OA=0.9I) 
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figure  31.      Computational  Results  for  the  proposed  Model  (OA=0.9I) 
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Figure  32.      Computational  Results  for  the  proposed  Model  (OA=  0.92) 
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Figure  33.      Computational  Results  for  the  proposed  Model  (OA=0.92) 
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34.      Computational  Results  for  the  proposed  Model  (OA=  0.93) 
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